On the Evolutionary Status of Class I Stars and Herbig-Haro 
Energy Sources in Taurus- Auriga 1 

Russel J. White and Lynne A. Hillenbrand 

Department of Astronomy, California Institute of Technology, MS 105-24, Pasadena, CA 

91125 

ABSTRACT 

We present high resolution (R ~ 34,000) optical (6330 - 8750 A) spectra 
obtained with the HIRES spectrograph on the W. M. Keck I telescope of stars 
in Taurus-Auriga whose circumstellar environment suggests that they are less 
evolved than optically revealed T Tauri stars. Many of the stars are seen only 
via scattered light. The sample includes 15 Class I stars and all Class II stars 
that power Herbig-Haro flows in this region. For 28 of the 36 stars observed, 
our measurements are the first high dispersion optical spectra ever obtained. 
Photospheric features are observed in all stars with detected continuum, 11/15 
Class I stars (42% of known Taurus Class I stars) and 21/21 Class II stars; strong 
emission lines (e.g. Ha) are detected in the spectra of all stars. These spectra, 
in combination with previous measurements, are used to search for differences 
between stars which power Herbig-Haro flows and stars which do not, and to 
reassess the evolutionary state of so-called protostars (Class I stars) relative to 
optically revealed T Tauri stars (Class II stars). 

The stellar mass distribution of Class I stars is similar to that of Class II 
stars and includes 3 spectroscopically confirmed Class I brown dwarfs. Class I 
stars (and brown dwarfs) in Taurus are slowly rotating (t>sim < 35 km/s); the 
angular momentum of a young star appears to dissipate prior to the optically 
revealed T Tauri phase. The amount of optical veiling and the inferred mass 
accretion rates of Class I stars are surprisingly indistinguishable from Class II 
stars. Class I stars do not have accretion dominated luminosities; the accretion 
luminosity accounts for ~ 25% of the bolometric luminosity. The median mass 
accretion rate of Class I and Class II stars of K7-M1 spectral type is 4 x 10 -8 
M Q /yr and the median mass outflow rate is 5% of the mass accretion rate. The 
large ranges in mass accretion rate (~ 2 orders of magnitude), mass outflow rate 
(~ 3 orders of magnitude) and ratio of these quantities (~ 2 orders of magnitude) 
represent real dispersions in young accreting stars of similar mass. We confirm 
previous results that find larger forbidden-line emission associated with Class I 
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stars than Class II stars. We suggest that this is caused by an orientation bias 
that allows a more direct view of the somewhat extended forbidden emission line 
regions than of the obscured stellar photospheres, rather than because of larger 
mass outflow rates. Overall, the similar masses, luminosities, rotation rates, mass 
accretion rates, mass outflow rates, and millimeter flux densities of Class I stars 
and Class II stars are best explained by a scenario in which most Class I stars 
are no longer in the main accretion phase and are much older than traditionally 
assumed. Similarly, although stars which power Herbig-Haro flows appear to have 
larger mass outflow rates, their stellar and circumstellar properties are generally 
indistinguishable from those of similar mass stars that do not power these flows. 

Subject headings: stars: formation — stars: low mass, brown dwarfs — stars: 
fundamental parameters — stars: winds, outflows — circumstellar matter - 
accretion, accretion disks 



1. Introduction 

Optical spectroscopy both initiated and continues to drive theories of low mass star 
formation. Strong HI and Call emission superimposed upon a late-type stellar absorption 
spectrum is one of the original defining characteristics of the T Tauri variable star class (Joy 
1945, 1949). Subsequent observations of these emission features has fueled theories for active 
chromospheres, mass accretion from disks, and mass loss in winds and jets (e.g. Calvet et al. 
1984; Lynden-Bell & Pringle 1974; Decampli 1981). High resolution optical spectroscopy 
is still the most accurate tool for characterizing the stellar properties (e.g. T e jj, logg, 
v sim, [Fe/H]) and investigating mass accretion and stellar jet/wind processes of young stars 
(Johns-Krull et al. 1999; Stempels & Piskunov 2003; Hartigan et al. 1995; Muzerolle et al. 
1998). For many young accreting stars, high resolution spectroscopy is often the only tool 
capable of extracting the stellar photospheric properties from the optically veiled spectrum 
(e.g. White & Basri 2003). However, because of high extinction, optical observers have been 
inhibited from pushing to the earliest stages of low mass star formation. These objects, often 
called protostars or Class I stars (Lada 1987), consist of a young star embedded within, 



1 The data presented herein were obtained at the W.M. Keck Observatory, which is operated as a scientific 
partnership among the California Institute of Technology, the University of California and the National 
Aeronautics and Space Administration. The Observatory was made possible by the generous financial 
support of the W.M. Keck Foundation. 
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but perhaps beginning to emerge from the collapsing envelope of material from which it is 
forming. Infall models with infall rates of a fewxlCT 6 M /yr are consistent with spectral 
energy distributions (SEDs) and scattered light images of Class I stars (e.g. Kenyon et al. 
1993a,b; Whitney et al. 1997, assuming a T Tauri-like mass and luminosity). This is 
believed to be the main accretion phase when the majority of stellar mass is acquired. How 
and how much of the infalling envelope material is transferred onto the central star, and 
the properties of this central star are largely unknown. Measuring stellar properties and 
the processes at the star-disk interface generally requires spectroscopic observations at short 
wavelengths (< 2/iin) where these self-embedded stars are the faintest. As a consequence, the 
current populations of protostars are all assumed to be younger than T Tauri stars because 
of their less evolved circumstellar environment; there is currently no photospheric evidence 
that demonstrates that these obscured stars are younger than T Tauri stars. 

The available (but very limited) spectroscopic observations of some embedded stars 
suggest that they are T Tauri-like (i.e. Class II-like), but with heavily veiled spectra and 
strong emission lines, implying high mass accretion rates and powerful stellar jets (Casali & 
Eiroa 1996; Greene & Lada 1996; Kenyon et al. 1998). The onset of a stellar jet or wind is 
generally believed to be the mechanism which clears the surrounding envelope, revealing the 
central star (e.g. Shu et al. 1987). At least ~ 60% of all Class I objects are associated with 
Herbig-Haro (HH) objects (Gomez et al. 1997; Kenyon et al. 1998; Reipurth 1999), which 
are regions of shocked gas believed to be powered by an outflow or jet (Reipurth & Bally 
2001). In comparison, only about 10% of Class II stars appear to power optical outflows, 
though this sub-sample has many of the highest accretion rates among T Tauri stars (e.g. 
DR Tau, DG Tau; Gullbring et al. 2000). Class II HH energy sources are often associated 
with spatially extended nebulosity (e.g. HL Tau, FS Tau; Stapelfeldt et al. 1995; Krist 
et al. 1998). Their semblance to younger protostars suggests that these stars could be 
transitioning between the Class I and Class II stages. 

A handful of studies over the last 2 decades have demonstrated using low resolution 
optical spectra that photospheric features can be seen in the scattered light continua of 
some Class I stars and HH energy sources (Mundt et al. 1985; Cohen et al. 1986; Kenyon 
et al. 1998). Although the advances in infrared techniques over this same time period led to 
a widespread shift from optical to infrared observations of the youngest stars (e.g. Greene & 
Lada 2000), optical spectroscopy still retains 2 advantages over near-infrared spectroscopy 
for characterizing the stellar and accretion properties. First, unlike near- infrared light which 
can be dominated by thermal emission from warm dust and gas, optical light is dominated 
by emission from the photosphere and the high temperature accretion shocks. It therefore 
offers the most direct view of the stellar properties and the accretion luminosity. Second, 
for small dust grains (< 1/iin), optical light scatters more efficiently than near-infrared 
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light. Thus, even if we can not see the stellar photosphere directly because of very high 
circumstellar extinction, the cavities commonly seen in the envelopes surrounding Class I 
stars (e.g. Padgett et al. 1999) may allow us to observe the photosphere and inner accretion 
processes through scattered light. This of course is only feasible in low column density star 
forming environments like Taurus-Auriga where the young stars are not deeply embedded 
within the large-scale molecular cloud. The protostars in Taurus present an opportunity to 
observe at optical wavelengths young stars at an unprecedentedly early age - potentially up 
to a factor of 10 younger than more optically revealed T Tauri stars if the statistical ages of 
embedded stars (Myers et al. 1987; Kenyon et al. 1990) are correct. 

In this paper we present a high spectral resolution optical survey of roughly 1/2 of 
all objects classified as Class I stars (the remaining 1/2 are currently inaccessible due to 
faintness) and nearly all objects (including both Class I and Class II stars) powering HH 
flows in Taurus. Our goal is to understand how stars evolve from embedded protostars to the 
optically revealed T Tauri stars. This requires determining if Class I stars are in fact younger 
than T Tauri stars, and if they are in the main accretion phase. Our approach is one which 
we hope will begin to bridge the observational segregation of T Tauri star and protostellar 
studies; although both the stellar and circumstellar properties of T Tauri stars are regularly 
investigated, the overwhelming majority of protostellar observations are directed towards 
understanding the circumstellar properties without regard for their stellar properties. 



2. Sample Definition, Observations and Data Analysis 

2.1. Sample 

The observational goal of this project was to observe a sample of stars in Taurus that are 
younger than T Tauri stars. Since the ages of young stars are often difficult to determine, 
especially at the earliest observable stages, we relied on 2 phenomenological indicators of 
youth to construct an initial target list. The first was based on measurements of the spectral 
energy distribution. We identified all stars that are Class I-like based on either their mid- 
infrared spectral index (Myers et al. 1987) or their bolometric temperature (Myers & 
Ladd 1993). The Infrared Astronomy Satellite (IRAS) provided an unbiased and complete 
(but flux and confusion limited) survey for such objects. The spectral indices, defined as 
a = dlog(XFx) / (d\og\) , were determined over the wavelength interval 2 to 25 /im using the 
K-[25] colors listed in Kenyon & Hartmann (1995). We included stars with a > 0, the 
generally accepted value used to distinguish Class I stars from Class II stars (Kenyon & 
Hartmann 1995). The bolometric temperatures, defined as the temperature of a blackbody 
with the same mean frequency as the observed SED, were taken from the compilation of 
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Chen et al. (1995). We included stars with T bo i < 650 K, the value proposed by Chen et 
al. (1995) to distinguish Class I stars from Class II stars. This sample consists of 32 stars. 
We caution that the spectral indices and the bolometric temperatures are in approximately 
half of the cases (57%) determined from SEDs that include flux contributions from multiple 
stars. As an extreme example, HH30 IRS, HL Tau, XZ Tau and LkHa 358 only have a 
combined, spatially unresolved 25 /zm measurement. In these cases the evolutionary classes 
of the components are assigned that of the integrated system. We note that if the 25 /im 
flux is equally divided among the components and individual indices are calculated using 
spatially resolved K-band magnitudes, these indices generally agree with the system's index 
(typically ±0.2 dex). However, there are remarkable exceptions. If the 25 /im flux of CoKu 
Tau 1 is divided in this way from its optically bright neighbor Hubble 4 (a weak-lined T 
Tauri star), its spectral index increase from —0.48 to +0.81, suggesting that it is a Class I 
star instead of a Class II star. 

The second indicator of youth we used is the presence of a spatially resolved optical jet. 
As discussed in the Introduction, stars which power HH flows are almost all either Class I 
stars or nebulous Class II stars. Adopting this selection criterion helped to identify a less 
SED-dependent sample; several apparently very young stars (based on spatially extended 
circumstellar emission) have insufficient photometric measurements to construct a SED (e.g. 
CoKu Tau 1, HV Tau C), and would have been excluded otherwise. Specifically, we included 
all stars that that are "Suspected Sources" of HH jets/flows in Taurus-Auriga as listed in 
the General Catalogue of HH Objects (Reipurth 1999). We also included the young star 
ZZ Tau IRS, which is the suspected energy source of HH 393 (Gomez et al. 1997), and has 
never been observed spectroscopically. We refer to stars that power a HH flow as HH stars, 
and those which do not as non-HH stars. This second selection criterion added 12 additional 
stars, for a total sample of 44 stars. 

In light of the circumstellar selection criteria used to identify this sample, we will refer 
to stars in this sample as environmentally young. We distinguish a strict sub-sample of 
bona fide Class I stars, by requiring that both a > 0.0 and 70 < T^ < 650° K (if only 
one of these SED measurements is available, the Class classification is based solely on that 
quantity). Stars with T bo i < 70° K are considered Class stars. Following this criterion, 
26 of the 44 environmentally young stars are Class I stars, and 1 is a Class star (IRAS 
04368+2557). The remaining 17 stars are classified as Class II stars. 

This combined sample of 44 stars is listed in Table 1 and is presumed to include the 
youngest stars in Taurus; the list does not include stars not detected at 2 /im (e.g. IRAM 
04191+1522; starless cores Motte & Andre 2001; Onishi et al. 2002). The components of 
binary systems (e.g. Ghez et al. 1993; Leinert et al. 1993; Duchene et al. 2004) are indicated 
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with the suffix A or B while the spatially unresolved multiples are indicated with the suffix 
AB or ABC. Stars newly observed in this study are distinguished (new) in the third column 
from stars for which we use previous high dispersion measurements (prev). Coordinates 
from the 2MASS database are given and in most cases are a substantial improvement to 
the available astrometry, especially for stars within dark cloud cores. We also list the 2 to 
25 /im spectral indices, bolometric temperatures from Chen et al. (1995), our evolutionary 
classification (I or II), whether the source powers a spatially resolved optical jet (Gomez et 
al. 1997; Reipurth 1999), if the star is visible in the POSS-II Red plates (inspection by eye), 
I c magnitudes and references, and 2MASS J-band magnitudes. 31 of the 44 stars power an 
optical jet; 15 of the 27 Class I/O stars have optical jets. 28 of these 44 stars are visible on 
the red POSS-II plates, and thus emit some observable light short-ward of 1 /im; the 16 stars 
not visible consist of 15 Class I stars and 1 (borderline) Class II star, IRAS 04154+2823. 

Figure 1 shows the distributions of a and T feo ; for the stars spectroscopically observed 
here (Section 2.3), and for the larger sample of Taurus T Tauri stars in Kenyon & Hartmann 
(1995) which have both measurements available. For illustrative purposes, the 3 Class I 
stars without T bol estimates (IRAS 04260+2642, IRAS 04325+2402, Haro 6-28) are plotted 
assuming a value of 500° K, and the 1 Class II star without a T bo i estimate (CoKu Tau 1) is 
plotted assuming a value of 2000° K. 

2.2. 7 c -band Imaging, Reduction, and Photometry 

Using the W. M. Keck II telescope and the Low Resolution Imaging Spectrograph (LRIS; 
Oke et al. 1995) in its imaging mode, we obtained 7 c -band images of 17 of the 44 environ- 
mentally young stars (mostly Class I stars) listed in Table 1 on 1998 Oct 30-31 and 1999 Dec 
13. In all cases, the integration time was set to 300 seconds which yields a sky count level 
of > 1/2 full- well. The images are 6x8 arcmin 2 in size and were processed by subtract- 
ing a median-averaged set of bias frames and dividing by a median-filtered set of flat field 
frames constructed from observation of the twilight sky. Equatorial "selected area" stan- 
dards (Landolt 1992) were observed for photometric calibration, assuming typical Mauna 
Kea extinction; all images were obtained at an airmass of < 1.1. Magnitudes were extracted 
using an aperture 6'/3 in diameter and a sky annulus extending from 26-30 arcseconds. These 
values are listed in Table 1. Errors are <0.1 magnitude unless marked with a colon. The 
measurements errors for stars listed with a colon are larger and more difficult to quantify 
because their emission is primarily spatially extended. 
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2.3. High Resolution Optical Spectroscopy 

The W. M. Keck I telescope and High-Resolution Echelle Spectrometer (HIRES; Vogt 
et al. 1994) were used by us to conduct observations on 3 separate observing runs: 1999 
Dec 5-6, 2002 Dec 13, and 2003 Feb 17-18. Variable and sometimes heavy cirrus conditions 
dominated the 1999 Dec run, while the other 2 runs had apparently photometric conditions. 
HIRES was used with the red collimator and the RG-610 filter. The Dl decker (1.15" x 14") 
was used, which has a projected slit-width of 4 pixels yielding a resolving power of ~ 34,000. 
The cross-disperser and echelle angles were set (at approximately 1.41 deg and -0.28 deg, 
respectively) to achieve a wavelength coverage of ~ 6330 to 8750 A, spanning 16 orders with 
wavelength gaps between the orders. This wavelength range includes several regions with 
temperature sensitive photospheric features and many of the strongest permitted emission 
lines (Ha, OI 8446, Ca II 8498, Ca II 8662) and forbidden emission lines ([OI] 6364, [Nil] 
6583, [SII] 6716, [SII] 6731) seen in T Tauri star spectra. Calibration frames obtained 
included those of an internal Quartz lamp for flat-fielding and a ThAr lamp for wavelength 
calibration. 

From the environmentally young star sample listed in Table 1, observational priority 
was given to targets which showed optical emission in the POSS-II Red plates (28/44) and 
for which no previous high resolution spectra were available. We observed 23 of the 28 
optically visible stars. The remaining five stars have been observed previously with either 
high dispersion optical spectroscopy (CW Tau, GK Tau, DO Tau, UY Aur; Hartigan et al. 
1995) or spatially resolved medium dispersion optical spectroscopy (XZ Tau A&B; Hartigan 
& Kenyon 2003). We include the available results of these previous measurements in our 
analysis and discussion (Sections 3 and 4). Finally, we observed 3 of the 16 environmentally 
young stars that are not visible on POSS-II plates, but are seen in our LRIS images (IRAS 
04154+2823, IRAS 04295+2251, IRAS 04489+3042), yielding a total of 26 of the 44 stars 
listed in Table 1, 31 including previous spectroscopy. 

In addition to these primary targets, we obtained spectra of 6 additional Taurus T Tauri 
star systems. IRAS 04278+2253 A&B, Haro 6-13, IRAS 04303+2240 all have bolometric 
temperatures less than 1000° K and have never been observed at high dispersion. HK Tau 
A&B have been classified as transitional Class I/II stars (Kenyon & Hartmann 1995); HK 
Tau B is also an edge-on disk system (Stapelfeldt et al. 1998; Koresko 1998) and has never 
been observed optically at high dispersion. We observed this pair of stars, even though the 
system's a and T& ; do not corroborate the claim that it is a transitional Class I/II star 
(a = —0.51, Tb i = 2148° K). The lightly veiled T Tauri stars DN Tau and V836 Tau were 
observed as well. The properties of these 6 star systems are also listed in Table 1, separately. 
They are included in the analysis, but distinguished from the environmentally young star 
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sample. 

The spectroscopic observations spatially resolved the components of the binary systems 
MHO 1/2 IRAS 04278+2253, HK Tau, and RW Aur. For the HV Tau system, only a 
spectrum of HV Tau C was obtained. In the case of L1551 IRS 5, we obtained spectra of the 
point-like optically visible object at the position of IRS 5, at the position of the HH object 
~ 12 arcseconds east of IRS 5, and at a position centered on the brightest portion of the 
optical reflection nebula that is 2'.5 directly west of IRS 5 (see Stocke et al. 1988). 

For 28 of the 36 stars observed (in 32 systems), these measurements are the first high 
dispersion optical spectra obtained. Exposure times and epochs of observations are listed in 
Table 2. To assist in the spectroscopic analysis, we observed 4 slowly rotating weak-lined T 
Tauri stars in the TW Hydrae association (TWA 8 A & B, TWA 9 A & B) and numerous 
dwarf and giant spectral type standards, 9 of which have accurately known radial velocities 
(Nidever et al. 2002). The observed spectroscopic standards were chosen to span a range 
of G, K and M spectral types at both dwarf and giant surface gravities. 

2.4. Spectroscopic Reduction 

The HIRES data were reduced using the facility "makee" reduction script written by 
Tom Barlow. The reductions included bias subtraction, flat-fielding, spectral extraction, sky 
subtraction, wavelength calibration, and heliocentric radial velocity corrections. The pro- 
jected spectra from target stars have typical full-widths at half-maximum of 1-2 arcseconds, 
though 5 stars (IRAS 04016+2610, IRAS 04302+2247, L1551 IRS 5, DG Tau B, Haro 6-5 B) 
have much broader profiles because of spatially extended emission. In these broader 
extraction aperture was used to include this extended emission; the sky level was determined 
outside these apertures. 

The 2 binary stars RW Aur and IRAS 04278+2253 were observed with the components 
aligned along the slit. The spectra of the components were extracted with a choice of 
apertures that minimize contamination from the companion. In the case of RW Aur, in 
which the secondary is separated by only 1'.'4 from the primary and is nearly 10 times fainter, 
the extracted spectrum of the secondary star is still contaminated by that of the primary; 
the secondary's spectrum showed a broad, but very distinct low level component to the Ca 
II emission lines (A 8498 A, A 8662 A), as well as very weak OI circumstellar absorption, 
both of which we attribute to contamination RW Aur A. These artifact features disappear 
cleanly if a spectrum of RW Aur A, scaled to 20 ± 5% of RW Aur B's continuum level, is 
subtracted from the spectrum of RW Aur B. Our final reduced spectrum of RW Aur B has 
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this 20% subtraction. 

Figure 2 shows portions of the HIRES spectra for the 28 stars with a continuum signal- 
to- noise ratio (S/N) greater than 2. Table 2 lists the continuum S/N per pixel at 6700 A and 
8400 A for all observations. For the 6 stars with S/N < 2 (2 Class II stars: DG Tau B and 
IRAS 04154+2823, and 4 Class I stars: MHO 1, L1551 IRS 5, IRAS 04295+2251, and IRAS 
04302+2247) the continuum detected (if real) is insufficient to measure any photospheric 
lines. These 6 stars are excluded from our photospheric analyses. Li I (6708A) is detected 
in the spectra of 24 of the 28 stars with a genuine continuum. Ha emission is detected in all 
36 stars observed. 



3. Analysis 

Figure 1 illustrates how this survey has extended the sample of stars with spectroscopic 
detections of the photosphere well into the Class I regime. Of the 26 Class I stars in Tau- 
rus (Table 1), 11 now have spectroscopic measurements of their photospheres and 15 have 
optical emission-line measurements. In addition to the remaining 15 Class I stars with no 
spectroscopic detection of their photosphere, there are 7 Class II stars from Kenyon & Hart- 
mann (1995) with no spectroscopically measured photospheres: IRAS 04370+2559, IRAS 
04414+2506, IRAS 04301+2608, FZ Tau, FT Tau, IRAS 04200+2759, and T Tau South. 
Note that only 3 of these Class II stars have both a and T hoi measurements and thus are 
shown in Figure 1. 

Though not complete, the considerable expansion of the available spectroscopy towards 
environmentally younger stars enables for the first time an investigation of the stellar and 
circumstellar properties between evolutionary classes. We wish to accomplish this without 
relying too strongly on any one evolutionary diagnostics (T^, a, HH properties), which may 
be flawed and/or have unaccounted for biases. Therefore we take the approach of comparing 
all measured and inferred properties with both the bolometric temperatures and spectral 
indices (when available), and distinguishing HH stars from non-HH stars. Possible differences 
between the properties of Class I stars and Class II stars (defined by T bo i and a criteria; 
Section 2.1) and HH stars and non-HH stars are investigated by comparing median values 
and by conducting K-S tests on these distributions (Press et al. 1997). K-S probabilities 
are reported in most cases; we consider distributions to be different if the K-S probability is 
less than 0.056 (ie. 2a difference). 

To increase the significance of these comparisons we identify and include an additional 
sample of Class II T Tauri stars. This sample includes the 6 additional Taurus T Tauri 



- 10- 



st ars observed here (all optically veiled), and the 26 optically veiled (at ~ 6000 A) Taurus 
T Tauri stars from the high spectral resolution studies of Basri & Batalha (1990), Hartigan 
et al. (1995) and White & Basri (2003). Our comparisons thus include only stars that 
are either environmentally young (Table 1) or are classical T Tauri stars defined by their 
optically veiled spectra. No weak-lined T Tauri stars (Class III stars) are included. 

3.1. Measured Spectral Properties 

3.1.1. Radial and Rotational Velocities 

Radial and projected rotational velocities are determined via a cross-correlation analysis 
(e.g. Hartmann et al. 1986). First, we identify up to 10 spectral regions of length 20-30 A 
with sufficient S/N (> 2) and without prominent telluric absorption features (e.g. Tinney 
& Reid 1998), strong gravity sensitive features (K I 7665 & 7699 A, Na I 8183 & 8195 A, 
Ca II 8498 A), or emission lines. These regions are cross-correlated with at least 3 spectral 
standards of spectral type similar to the target star (see Section 3.1.2) and with radial 
velocities accurate to 0.3 - 0.4 km/s (Nidever et al. 2002). 

To determine the final radial velocities of the stars, the velocity offsets from the stan- 
dards must be corrected for barycentric velocity corrections and possible errors in the wave- 
length solution of either the standard or target spectrum. The former is determined using 
the rvcorrect task in IRAF. The latter correction may only be relevant when combining data 
from multiple nights and/or multiple runs, as is the case here. This correction is determined 
by cross-correlating the telluric A-band of the standard and the target star. Failure to cor- 
rect for this effect leads to radial velocity errors of 3-4 km/s. Correcting for this effect yields 
consistent radial velocities of the standards to better than 0.5 km/s, consistent with the 
accuracy of their radial velocities. Uncertainties are estimated by adding, in quadrature, the 
standard deviation of the multiple radial velocity estimates (3-7 standards are used in each 
case) and the uncertainty in each radial velocity estimate, which is assumed to be the error 
in the mean of all spectral regions used per star. 

Rotational velocities (vsmi) for the observed targets are determined from the width of 
the peak in their cross-correlation with a slowly rotating standard. The width is estimated 
by fitting a Gaussian plus quadratic function to the cross-correlation peak. The width 
versus -usim relation is determined empirically by cross-correlating a "spun-up" spectrum of 
a slowly rotating standard. The rotationally broadened spectra are constructed using the 
profiles given in Gray (1992; e = 0.6). As has been noted in previous work (e.g. White 
& Basri 2003), the presence of a continuum excess (see Section 3.1.3) has little effect 
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on the inferred projected rotational velocities, as long as the features are sufficiently well 
measured. Uncertainties in vsini are determined following the same procedure used for the 
radial velocities. We find that all templates within a few spectral subclasses of the target 
star generally give consistent vsini values; an exact match in spectral type is not critical 
in determining vsini. However, differences in spectral type and luminosity class can affect 
vsini measurements at low rotational velocities (less than ~ 10 km/s; see White & Basri 
2003). To avoid introducing this bias, we adopt the conservative approach in determining 
vsini upper limits. If the uncertainty is less than 1 standard deviation from the theoretical 
vsini detection limit (7 - 15 km/s, set by S/N), the measurement is assigned an upper limit. 

The described procedure worked well for determining radial and projected rotational 
velocities for 24 of the 28 stars with a continuum detection. Because of low S/N and, in 
some cases, high veiling (Section 3.1.2), the cross-correlation analysis for the remaining 4 
stars (IRAS 04016+2610, IRAS 04264+2433, IC 2087, and IRAS 04489+3042) was successful 
in only a few segments (2-3) of their spectra. The correlation results for these 4 stars imply 
no measurable rotation, but do yield radial velocities that are consistent with that of Taurus. 
In the cases of IRAS 04016+2610 and IC 2087, we are especially cautious about the results. 
The photospheric-like features in their spectra are quite weak and/or noisy, and possibly 
could be spurious noise events. The results for IRAS 04264+2433 and IRAS 04489+3042 
are slightly more robust, however, since both stars show TiO absorption bands that can be 
more confidently identified (see Figures 1 & 2). The radial velocities and vsini detection 
limits (15 km/s) for these stars are listed in Table 2 followed by colons to indicate possible 
larger systematic uncertainties. 

Although for the majority of these stars their measurements are the first epoch of 
high resolution spectroscopy, no definitive spectroscopic binaries are identified. With the 
exception of HV Tau C, all stars have radial velocities within 2a of the mean of Taurus- 
Auriga (17.4 ± 2.1; Hartmann et al. 1986), and thus are consistent with being Taurus 
members. The radial velocity of HV Tau C (23.1+0.6) is 2.6cr from the mean of Taurus, 
and thus the star could be a single-lined spectroscopic binary. Unfortunately we do not 
have a spectrum of its companions HV Tau AB for a more strict comparison. We also note 
that radial velocities of HK Tau A and HK Tau B are significantly different (21.6 ± 0.5 and 
17.4 ± 0.3), several times that expected from relative orbital motion alone. We suggest that 
either one of these stars is a single-lined spectroscopic binary or that the pair is not physically 
associated. The latter case would explain their apparently non-coplanar disks (Duchene et 
al. 2003). 

The projected rotational velocities of the environmentally young star sample are rela- 
tively slow, ranging from 10 - 34 km/s. The measured vsini values of this sample and the 
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additional sample of Class II stars are plotted in Figure 9 versus the evolutionary diagnostics 
T bo i and a. The environmentally young stars have measured t>sim values (excluding upper 
limits) that are generally quite similar to the values of environmentally older stars. Treating 
i>sim upper limits as detections, the distributions of vsmi values for Class I and Class II stars 
are similar (K-S probability of 0.079), though the median vsmi of Class I stars is slightly 
larger than that of Class II stars (median Class I fsim = 18.0 km/s, o = 5.1; median Class 
II vsmi = 14.4 km/s, a = 10.8). Comparing HH stars and non-HH stars, the distributions 
of fsim are more discrepant (K-S probability of 0.0034; a 2.9a difference); HH stars rotate 
slighly faster than non-HH stars (median vsmi = 19.4 km/s, a = 5.9; median vsmi = 12.1 
km/s, a = 11.5). Table 3 summarizes the median statistics, excluding brown dwarfs. 



3.1.2. Spectral Types and Continuum Excesses 

Spectral types and continuum excesses are determined by matching the line/band ratios 
and line/band depths, respectively, to those of rotationally broadened and artificially veiled 
dwarf spectral standards. The procedure works as follows for a given star. First, all of the 
comparison dwarf spectral standards are rotationally broadened to have the same t>sim as the 
young star. For each rotationally broadened standard, a set of "veiled" templates are then 
generated by adding a range of continuum excess values, all assumed to be constant over each 
~ 20 A spectral region. This set of rotationally broadened and veiled templates are then 
compared to the young star sample to determine which spectral type gives the best match. 
Quantitatively, the best match is determined by minimizing the root-mean-squared values 
over 3 temperature sensitive regions: 6490 - 6502 A, 7120 - 7128 A, and 8431 - 8436 A (see 
Figures 2-4). In most cases, the assigned spectral type is that of the spectral standard which 
gives the best fit. In some cases, however, emission lines or anomalous low S/N features yield 
clearly discrepant spectral types. Consequently, all spectral types are confirmed by visual 
inspection. Uncertainties in the spectral type are determined by the range of spectral types 
which can reasonably match the spectrum; the uncertainties range from 0.5 to 3 spectral 
subclasses. More robust techniques for determining temperatures such as synthetic modeling 
of the spectra are currently unwarranted given the low S/N of the data. The low S/N of 
IRAS 04016+2610 and the high optical veiling of IC 2087 prevent us from determining their 
spectral types. However, both show Ti I (8435 A, 8436 A) absorption that is stronger than 
Nb I (8439 A) absorption, suggesting a spectral type cooler than G. The absence of any clear 
TiO bands suggests a spectral type hotter than M. We classify these two stars as K stars. 
The inferred spectral types for all stars are listed in Table 1. 

The inferred spectral types range from K0 to M6 for the environmentally young star 
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sample. Class I stars and Class II stars have a similar distributions of stellar spectral types. 
HH stars, on the other hand, tend to be slightly hotter than non-HH stars. Of the 46 stars 
studied here with spectral types hotter than M3.5 (M > 0.24 M ), 20 are HH stars (43%). 
In contrast, of the 9 stars with spectral types of M3.5 or cooler, only 2 are HH stars (22%). 

Once the spectral type of the star is determined, the continuum excess, called r and 
defined as r = F excess / F photosphere , is determined by comparing the depths of photospheric 
features to those of a rotationally broadened and artificially veiled dwarf spectrum of the 
same spectral type. The continuum excess added to the standard is allowed to be negative; 
this results in a comparison standard with deeper absorption features. Consequently, stars 
with no continuum excess should have r's that scatter about zero, and this scatter can be 
used to estimate the uncertainty in r. The continuum excess is determined in 2 spectral 
region, 6500 A and 8400 A, by averaging the results within 50 A of these wavelengths. The 
uncertainties in the continuum excess are estimated by changing the rotationally broadened 
and veiled template by the spectral type uncertainty and redetermining the continuum excess; 
the spectral type uncertainty typically dominates the uncertainty in the continuum excess. 
If the uncertainty in the continuum excess is larger than the measured value, the uncertainty 
is assigned as an upper limit to detection. As a check of our procedure, we note that the 
4 weak-lined T Tauri stars TWA 8a, TWA 8b, TWA 9a, and TWA 9b, with spectral types 
ranging from K7 to M5.5, can be successfully modeled using a rotationally broadened dwarf 
spectrum without any continuum excess (e.g. < r 650 o > = —0.04 ±0.11). For the K stars 
IRAS 04016+2610 and IC 2087, we use a K4 dwarf as the template from which the veiling is 
determined. The inferred continuum excess values at 6500 A (r 650 o) and at 8400 A (r 840 o) 
are listed in Table 2. For 5 stars, the S/N is insufficient to determine the continuum excess 
at 6500 A but is sufficient to measure the excess at 8400 A. 

Figure 10 shows the ratio of the continuum veiling at 8400 A to the continuum veiling at 
6500 A as a function of stellar effective temperature. Values plotted include measurements 
from this study (Table 2) and from White & Basri (2003). In addition, we show the ratio of 
optical veiling values from Hartigan & Kenyon (2003), determined from similar but slightly 
shorter wavelength intervals (r6ioo/ r 8ii5), using low resolution spectra; only veiling values 
above 0.1 are included from this data set to avoid spuriously large ratios. For early- to 
mid-K stars, the continuum excesses at 8400 A are in many cases larger than the values 
determined at 6400 A. This is somewhat surprising, since the optical excesses are generally 
attributed to featureless continuum emitted by a ~ 10 4 K hot spot(s) (e.g. Kenyon et al. 
1994; Calvet & Gullbring 1998), and therefore should be significantly less at 8400 A. We 
note that this effect is consistent with previous veiling studies which find the veiling to be 
"almost flat" in the red (Basri & Batalha 1990); this previous work primarily focused on 
K5-M0 type stars for which an excess of constant flux would yield nearly a constant veiling 
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value (r 8 4oo/^650o ~ 1-0; Figure 10). We postpone discussion of this interesting phenomenon 
until Section 3.3.1, but use these empirical ratios to estimate the continuum excess at 6500 
A for the 5 stars with a r 840 o measurement, but no r 650 o measurement. For the 2 K stars, 
r65oo is estimated to be rs4oo/l-2; for the Ml star, ^500 is estimated to be to be rg4oo/0.70; 
for the M5.5 and M6 stars, r 650 o is estimated to be r 8 4oo/0.43. These estimated r 650 o values 
are listed in Table 2 followed by colons. 

For IRAS 04303+2240, which was observed both on 1999 Dec 6 and 2003 Feb 17, no 
photospheric features were detected during the first epoch, but early-M spectral features 
were clearly identified during the second (Figure 1). The star was also much brighter during 
the first epoch than the second epoch; the relative counts in the spectra suggest a brightness 
decrease by factors of 3.25 at 6500 A and 2.47 at 8400 A (see S/N in Table 2). Although 
no continuum excesses could be measured in the first epoch since the photosphere is not de- 
tected, we estimate the excesses during this epoch by assuming that the change in brightness 
between the two observations is caused only by a change in the continuum excess. These 
larger first epoch values are listed in Table 2 and followed by colons. 

The distribution of continuum excesses at 6500 A (r 650 o) is shown in Figure 11 as a 
function of T bo i and a. Values for the additional sample of Class II stars are also shown, 
with continuum excess measurements from Basri & Batalha (1990); Hartigan et al. (1995) 
and White & Basri (2003). Environmentally young stars have optical veiling values that 
are similar to those of environmentally older stars. K-S tests show no differences in the 
distributions, but slight differences are seen in the median values. The continuum excesses 
of the 11 Class I stars (median rgsoo — 1-0, cr = 0.6) are marginally larger than the continuum 
excesses of the 50 Class II stars (median r 650 o = 0.6, a = 1.3). Similarly, the continuum 
excesses of the 24 HH stars (median r 650 o = 11, o = 1.3) are marginally larger than the 
continuum excesses of the 36 non-HH stars (median r 650 o = 0.6, cr = 1.0). Table 3 summarizes 
the median statistics, excluding brown dwarfs. 



3.1.3. Surface Gravity Indicators 

Pressure sensitive photospheric features offer a means of determining if the surface 
gravities of Class I stars are less than those of Class II stars, for a given mass. This would 
indicate that they have larger radii and hence younger ages. Some of the more distinguishing 
pressure sensitive features for K and M spectral types include CaH (6382 A, 6389 A, 6903 - 
6946 A), K I (7665 A, 7699 A), Na I (8183 A, 8195 A), Ca II (8498 A, 8542 A, 8662 A) and 
Fe I lines (see e.g. Kirkpatrick et al. 1991). The molecular absorption bands of TiO and VO 
present in M dwarfs are also known to be pressure sensitive (Torres-Dodgen & Weaver 1993; 
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White & Basri 2003). Our spectra include K I (7699 A), Ca II (8498 A, 8662 A), several 
Fe I lines and CaH and TiO absorption bands. Of these, Ca II is not a useful diagnostic of 
surface gravity, since it often shows either strong emission from accretion activity, or core 
emission from chromospheric activity. 

The spectra shown in Figures 2-4 include the K I at 7699 A. This feature is one of the 
strongest and most easily identifiable absorption features in our modest S/N optically veiled 
spectra. However, many stars show non-photospheric components to this feature, including 
both broad emission profiles with narrow absorption superimposed (e.g. DG Tau, HL Tau, 
T Tau) and broad asymmetric absorption with possible emission (e.g. IRAS 04303+2240, 
IC 2087, RW Aur A). For objects with no obvious non-photospheric emission, however, the 
wing strengths of K I for both Class I and Class II stars are intermediate between those of 
narrow giants and those of pressure broadened dwarfs. This is a well known characteristic 
of T Tauri stars. 

For early- and mid-M spectral types, we find that dwarfs, as opposed to giants or some 
combination of the two, provide a better match to the pressure sensitive Fe I lines and TiO 
bands of both Class I and Class II stars White & Basri (see also 2003). In addition, CaH 
bands at 6382 A and 6389 A are clearly seen in many of the Class I and Class II stars 
observed here. These features are strongest in the spectra of early M dwarfs, and absent in 
spectra of similar spectral type giants. 

The combined results of these comparisons suggest that the surface gravities of both 
Class I and Class II stars are intermediate between those of dwarfs and giants, but somewhat 
more similar to dwarfs. This is consistent with the predicted surface gravities for young stars 
(e.g. Siess et al. 2000), whose radii should be only a few times their main sequence values, 
as opposed to ~ 100 times that of their main-sequence values in the case of giant stars. 
Theoretical predictions suggest Class I stars should have radii less than 2 times those of T 
Tauri stars (e.g. Stahler 1988), which translates into differences in surface gravity of less than 
a few tenths of a dex. Given the relatively low S/N of the Class I spectra, measurements of 
these small log(g) differences are not possible from line profile analysis. We are also cautious 
of such an approach, given the non-photospheric emission associated with gravity sensitive 
features (K I 7699 A; Section 3.1.4). Analyses which focus on determining surface gravities 
from molecular-band strengths (e.g. Mohanty et al. 2004) offer a more promising approach. 
With the current data at hand, we can neither prove nor disprove that Class I stars have 
lower surface gravities, as would be expected if they are younger than Class II stars. 
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3.1.4- Prominent Emission Lines 

The emission-lines of young stars are powerful probes of the mass accretion and mass 
outflow processes. Their ratios characterize temperatures and densities of the emitting re- 
gions and their intensities and velocity profiles trace the flow of circumstellar material (e.g. 
Edwards et al. 1987; Hamann & Persson 1992; Hamann 1994; Hartigan et al. 1995). In 
Figures 5-8 are shown emission-line profiles of the observed Taurus stars in Table 2. The 
plotted profiles include the 2 permitted lines Ha 6563 A and Call 8498 A, and the 3 for- 
bidden emission-lines [01] 6364 A, [Nil] 6583 A, and [SII] 6731 A. Fe II 6456 A, K I 7699 A 
and 01 8446 A emission can be seen in the spectra plotted in Figures 2-4, when present. 

Equivalent widths (EWs) of Ha, Call 8498 A, Call 8662 A, [01] 6364 A, [Nil] 6583 
A, [SII] 6716 A, and [SII] 6731 A are listed in Table 3. These EWs include emission at all 
velocities which, for several forbidden lines, means summing distinct velocity components. In 
cases where the Call emission appears in the center of the photospheric absorption feature 
(CoKu Tau 1, Haro 6-28, Haro 6-33, RW Aur B, HK Tau A & B, DN Tau, V836 Tau), 
the EWs are estimated from the approximate bottom of the stellar absorption profile; no 
corrections for the stellar absorption are made for stronger emission-line stars. Detection 
limits are estimated based on the S/N of the spectra and/or the quality of the sky subtraction 
for telluric sensitive features (e.g. [01] 6364 A). In cases where an emission-line is present, 
but the S/N in the continuum is less than 1.0, the equivalent width is assumed to be a 
lower limit (up to 999 A). Lower limit EW measurements of [SII] 6716 A and 6731 A use 
the same continuum value to preserve the ratio of these neighboring features. For stars 
with S/N < 1.0 and no emission-lines, no limits are given. Uncertainties in the EWs are 
generally 10% of the listed value, and usually limited by the difficulty in determining the 
continuum level of a cool type photosphere underneath an emission-line with broad wings. 
The EWs of [OI] 6364 A and to a lesser extent [SII] 6716 A have larger uncertainties (10- 
25%) because of the uncertain removal of sky emission and stellar absorption features at or 
near those wavelengths. In several cases, poor sky subtraction yielded lower limits on the 
EW measurements. 

Ha Emission 

Broad Ha emission is one of the most ubiquitous characteristics of mass accretion. Mea- 
surements of the line profiles potentially offer a means to characterize the accretion velocity 
and mass accretion rate (e.g. Muzerolle et al. 1998). Ha emission is detected from all stars 
observed, albeit marginally for IRAS 04154+2823 and IRAS 04295+2251. Figure 12 shows 
the distributions of Ha EWs versus T bo i and a for the environmentally young star sample 
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and the additional Class II star sample. For illustrative clarity, no lower limits are plotted 
(typical of stars with no continuum detection); since all stars show Ha emission, there are 
no upper limits. Environmentally young and old stars have indistinguishable Ha EW dis- 
tributions, according to K-S tests. The median EW[Ha]s of the 16 Class I stars and the 49 
Class II stars are both —59 A. The median EW[Ha]s of the 25 HH stars is —63 A and is 
similar to that of the 39 non- HH stars, -58 A. Table 3 summarizes the median statistics, 
excluding brown dwarfs. 

Several environmentally young stars show evidence of a blue-shifted absorption compo- 
nent superimposed upon the Ha emission profile. This feature is often seen in the profiles of 
accreting TTSs and is usually attributed to absorption from an outflowing jet or wind (e.g. 
Edwards et al. 1987; Alencar & Basri 2000). 

The full widths at 10% of the peak flux levels (called 10%-widths, hereafter) of the 
Ha profiles are listed in Table 4. This profile measurement has been proposed as a useful 
diagnostic of accretion (e.g. White & Basri 2003). For the environmentally young star 
sample, these values range from 139 km/s to 610 km/s. The majority of the Ha profiles are 
broad (25/28 have 10%-widths > 250 km/s), consistent with the predictions of a magnetically 
channeled accretion flow (e.g. Hartmann et al. 1994). Nevertheless, several optically veiled 
and presumably accreting stars show relatively narrow 10%-widths (e.g. IRAS 04260+2642 
and HK Tau B). The distribution of Ha 10%-widths of Class I stars are marginally different 
for that of non-HH stars (K-S probability = 0.062; a 1.9c difference); Class I stars have 
narrower 10%-widths, in the median, than Class II star (356 km/s versus 429 km/s). In 
contrast, distributions of 10%-widths for HH stars and non-HH stars are indistinguishable 
(K-S probability = 0.77), and the median values are essentially the same (410 km/s and 413 
km/s, respectively). Table 3 lists the median values, excluding brown dwarfs. 

Forbidden-line Emission 

The optically thin forbidden emission-lines are believed to originate in an outflowing jet or 
wind. Their intensity is therefore expected to be directly proportional to the amount of 
material being funneled along the jet, as viewed through the slit of the spectrograph. The 
strongest forbidden emission-lines in our spectra are [OI] 6364 A and [SB] 6731 A, both of 
which are detected in 21 of the 28 environmentally young stars observed. Figure 13 shows 
the distributions of EW[SII]s versus T hoi and a. No lower limits are plotted, but upper limits 
are. Additional EW measurements of Taurus T Tauri stars are included using the values 
measured by Hartigan et al. (1995, all velocities). 

The environmentally young stars have EW[SII] within the range of values measured 
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for environmentally older stars, but the distributions are different. The median EW[SII] of 
Class I stars is larger than that of Class II stars (median Class I EW = —3.1 A; median 
Class II EW = —0.15 A) and the distributions are statistically different (K-S probability 
of 1.9 x 10~ 5 ). Similarly, the median EW[SII] of HH stars is larger than that of non-HH 
stars (median HH EW = —1.8 A; median non-HH EW = —0.14 A) and the distributions 
are different, though less significantly (K-S probability of 0.0086). Table 3 lists the median 
values, excluding brown dwarfs. These comparisons may be somewhat biased by the large 
fraction of [SII] non-detections (31% for Class I stars, 34% for Class II stars, 15% for HH 
stars, 47% for non-HH stars), which are treated as detections. Nevertheless, comparing 
only systems with detected emission, the EW[SII]s of environmentally young stars are still 
systematically larger than those of environmentally older T Tauri stars. 

The forbidden emission-lines of T Tauri stars are often blue-shifted and occasionally 
show double-peaked profiles, with one peak near the stellar velocity and the second offset 
blue- ward by 50 - 300 km/s (e.g. Appenzeller et al. 1984; Mundt et al. 1987; Edwards 
et al. 1987). Only 3 of the 21 environmentally young stars with forbidden-line emission 
shows a distinct double peaked emission profile (DG Tau B, L1551 IRS 5, HL Tau). One 
of the additional TTS stars observed also shows a double-peak (Haro 6-13). In the 2 cases 
for which the systemic velocity is known (HL Tau and Haro 6-13), the blue-shifted peaks 
are offset by approximately -200 km/s and -100 km/s, respectively. About a third of the 
21 environmentally young stars with [OI] emission have forbidden emission peaks close to 
the systemic velocity, but with wings extending blue-ward by > 50 km/s. Half of the stars 
observed, however, show only forbidden emission at or very near the systemic velocity. This 
fraction is slightly higher if only Class I stars are considered. Of the 13 Class I stars with 
[OI] emission, only 5 show evidence of a blue-shifted component (MHO-2, GV Tau A, L1551 
IRS 5, HL Tau, IRAS 04489+3042). This effect could be explained if the jet emission in the 
majority of Class I stars observed propigates perpendicular to our line of sight. 

The Effect of Orientation on the Emission-line Profiles 

The orientation of each star/disk/jet system relative to our line-of-sight may have effects on 
both the measured EWs and emission-line profiles. One of the most useful qualities of an 
EW is that, by definition, it offers a direct measure of the emission-line flux relative to the 
continuum, independent of extinction. Thus, if the spectra are not flux calibrated and the 
extinction is not accurately known, as is the typical case here, the total emission-line flux may 
still be determined if the intrinsic flux from the stellar continuum can be estimated (e.g. by 
assuming a flux density distribution based on age and distance) . However, this assumption is 
only valid if the continuum and emission-line fluxes originate from the same region, and thus 
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have the same extinction. This is most likely true for the permitted emission-lines of young 
stars, which in general appear to be optically thick and originate in high density regions 
close to the star. However, this assumption may not be true for the forbidden emission lines, 
which often originate in regions that are spatially extended from the star (e.g. Reipurth 
1999). In such a case, the emission-line region may be more directly observable than the 
young partially embedded central star is. The preferentially attenuated continuum flux will 
consequently produce artificially large EW values. 

The three systems with well defined edge-on disks based on direct imaging, HH 30, HV 
Tau C, and HK Tau B, likely suffer from this bias to some degree. The forbidden emission 
line profiles of the edge-on disk systems are narrow and centered on the systemic velocity, 
consistent with a bipolar flow perpendicular to our line-of-sight. High spatial resolution 
Hubble Space Telescope images of both HH 30 and HV Tau C show spatially extended [01] 
and [SII] emission (Burrows et al. 1996; Bacciotti et al. 1999; Stapelfeldt et al. 2003); this 
emission is more directly observable than the continuum stellar emission which is scattered 
into our line of sight. Preferential attenuation of the continuum likely explains why the EWs 
of [SII] 6731 A for HH 30 and HV Tau C are more that 2 orders of magnitude larger than those 
of any Taurus T Tauri star in the survey of Hartigan et al. (1995). The possibly biased EWs 
of edge-on disk systems are therefore distinguished in Figure 13. A preferentially attenuated 
continuum may also explain the large forbidden emission-line EW lower limits of the 2 stars 
with no detected continuum (DG Tau B and L1551 IRS5). Without spatial information on 
the origin of these emission-line features and the efficiency of scattered continuum emission, 
this is not something we can correct for. We highlight this effect, however, as a potential bias 
of the forbidden emission-line EWs and the mass outflow rates inferred from them (Section 
3.3.2). 

In contrast to the forbidden emission-line EWs, the Ha EWs of the 3 edge-on type disk 
systems are consistent with the EW distribution of the larger, presumably randomly oriented 
sample. This agrees with the expectation that the majority of the Ha emission originates very 
close to the star and therefore experiences the same attenuation as the stellar continuum. 
However, they nevertheless have relatively narrow Ha 10%-widths (293 km/s, 261 km/s, 
and 196 km/s). On average these values are more consistent with those of non-accreting T 
Tauri stars, with typical 10%-widths < 270 km/s (White & Basri 2003). Assuming that 
the accretion occurs via the same magnetically chanelled high velocity flow for these 3 stars, 
then either the edge-on orientation prevents a direct view of the high velocity emission or 
the projected velocities along our line of sight are low. The absence of strong blue-shifted 
absorption component superimposed on the Ha emission is also consistent with an outflow 
perpendicular to our line of sight. 
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In light of the emission-line properties of the known edge-on disk systems HH 30, HV 
Tau C, and HK Tau B, we suggest that 3 additional stars with similar properties, CoKu 
Tau 1, IRAS 04260+2642, and ZZ Tau IRS, may also have nearly edge-on orientations. This 
prediction is based on (1) the relatively narrow Ha 10%-widths (285 km/s, 139 km/s, 268 
km/s), (2) the absence of any superimposed absorption on the Ha emission profiles, and (3) 
the unusually large strengths of their forbidden emission-lines, which are both narrow and 
very close to the systemic velocity (see Figures 5-7 and Table 3). The emission-lines of IRAS 
04158+2805 and IRAS 04248+2612 are also marginal consistent with an edge-on orientation, 
but the narrow profiles of these two stars may simply result from the lower infall/outflow 
velocities expected in substellar objects (e.g. Jayawardhana et al. 2003). In a few of these 
cases, however, we note that the narrow forbidden emission lines are centered ~ 10 — 40 
km/s blueward of the systemic velocity (e.g. IRAS 04158+2805, IRAS 04260+2642, ZZ Tau 
IRS), suggesting that these systems are close to, but not completely edge-on. 

Emission-line Variability 

L1551 IRS 5 and IRAS 04303+2240 show significant changes in emission-line properties over 
a time scale of ~ 3 years. During the first epoch observations of IRAS 04303+2240, when the 
star is much more optically veiled and presumably accreting at a higher rate, the permitted 
emission-lines (Ha, Ca II) show evidence of absorption from a blue-shifted outflow, but the 
forbidden emission-lines are very weak or absent. During the second epoch, when the star is 
less optically veiled, the permitted line profiles are slightly narrower, more symmetric, and 
show no evidence a superimposed absorption. The forbidden emission-lines, however, are 
considerably stronger. Even after accounting for the diminished continuum excess, the EW 
fluxes of [OI], [Nil] and [SII] are all an order of magnitude larger during the second epoch. 

No continuum emission is detected in the first observation of L1551 IRS 5, but is 
marginally detected in the second (S/N = 1.4 at 8450 A). Like IRAS 04303+2240, the 
emission-lines profiles of LI 551 IRS 5 changed significantly over the same 3- year time-scale. 
The changes included broader emission profiles for both Ha and the forbidden emission 
lines, the disappearance of the Ca II emission, and possibly weaker line-emission overall 
(lower upper-limits). It is interesting that the emission-line profiles of the first epoch of IRS 
5 are quite similar to the emission-line profiles of the HH knot obtained during the second 
epoch. The similarities include both the approximate velocity profile of the Ha emission and 
the 3-peaked profiles of the forbidden-line emission. 

The spectrum of the L1551 nebula is the only spectrum from this system that shows 
Ha emission with a peak close to the expected systemic velocity (i.e. the mean of Taurus). 
This may actually be scattered emission from the star. Low resolution spectra of this nebula 
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have identified stellar absorption features of a G-type photosphere (Mundt et al. 1985). No 
photospheric features were detected in our spectrum. Observations of this scattered light 
nebula may be the only way to observed the stellar emission from L1551 IRS 5 at optical 
wavelengths. 

3.2. Inferred Stellar Properties 

3.2.1. Stellar and Bolometric Luminosities 

One of the main motivations for determining the stellar luminosities of the stars identi- 
fied as environmentally young in Table 1 is to see if they are in fact younger than more op- 
tically revealed T Tauri stars. Since low mass stars contract primarily along Hayashi tracks, 
maintaining roughly constant temperature, younger stars should have larger radii and be 
more luminous. The luminosity of a young star is typically calculated from a reddening cor- 
rected photometric measurement (e.g. I c or J) and a bolometric correction corresponding to 
its spectral type. Although most of the stars observed here have both measured photometry 
and now a spectral type, ascertaining the intrinsic stellar energy distribution and the proper 
reddening correction necessary to estimate their stellar luminosities is inhibited by the exten- 
sive non-photospheric emission associated with these stars. At optical wavelengths, scattered 
light emission can increase the apparent stellar flux measured in a fixed aperture and can 
cause underestimates of extinction because of preferential blue scattering. At near-infrared 
wavelengths, the general trend of increasing thermal emission from the disk toward longer 
wavelengths can lead to over-estimates of the extinction. Despite these possible systematic 
uncertainties, we explore several traditional means for calculating stellar luminosities. We 
consider de-reddening the stars using colors which are as uncontaminated, that is as close to 
photospheric as possible, namely I c — J, J — H, and H — K s . I c magnitudes are listed in Table 
1; J, H, and K s magnitudes are from the 2MASS database. The underlying photospheric 
colors are adopted according to our spectral types and a standard interstellar reddening law 
(Cohen et al. 1981). For comparison purposes, we apply bolometric corrections to 3 de- 
reddened magnitudes: J c , J, and K s . For the fraction of our sample where all three methods 
can be applied, comparison of the results reveals systematic effects. Luminosities calculated 
from K s and H — K are systematically higher than those computed from J and J — H (by 
0.35 dex) or from I c and I c — J (0.45 dex). This is likely due to the systematically larger 
excess emission at K s caused by thermal emission from hot dust and gas. However, the 
differences in luminosities calculated from J, J — H and I c , I c — J are scattered about unity 
though with large dispersion. This suggests that either may provide an adequate measure 
of the stellar luminosities. Since we have J and H magnitudes for the majority of stars, we 
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adopt the stellar luminosities calculated from J — H extinctions and a bolometric correction 
to J-band. The calculated visual extinction and stellar luminosities are listed in the fourth 
and fifth columns of Table 5. The uncertainties in these estimates are dominated by system- 
atic effects, which lead to errors of 1-2 magnitudes in Ay and factors of 2 — 3 in L star , with 
even larger errors for edge-on disk systems (e.g. HH 30). 

The average visual extinction of the 11 Class I stars with stellar spectral types is 10.5 
mag, with a standard deviation 4.2 mag. This is only a factor of a few larger than typical 
values for Class II stars in Taurus, and is a factor of a few less than is often suggested by the 
rising energy distributions (Greene & Lada 1996) and scattered light images (e.g. Kenyon 
et al. 1993b; Whitney et al. 1997) of most Class I stars. As described in Section 2.1, our 
Class I sample is biased towards those which are more optically revealed. 

Figure 14 shows the distribution of stellar luminosities for environmentally young stars 
in Taurus as a function of effective temperature. The temperatures are estimated from 
the new spectral types listed in Table 2, assuming the dwarf temperature scale adopted by 
Hillenbrand & White (2004). The 6 environmentally young stars which are likely edge-on 
disk systems (Section 3.1.4) are marked since this orientation likely leads to underestimates 
of their stellar luminosity. The 3 panels separate stars with T bo i < 650, stars with a > 0.0, 
and stars which power HH flows. Also shown for comparison are the mean stellar luminosities 
of Class II and Class III T Tauri stars in Taurus (environmentally young stars are excluded) 
measured in 6 temperature bins. The T Tauri stellar luminosities are taken from Kenyon 
& Hartmann (1995), and are determined from a bolometric correction to the J-band, as 
was done here 2 . Figure 14 also shows the pre-main sequence evolutionary models of Siess 
et al. (2000). In comparison to these models, T Tauri stars in Taurus have a mean age 
of a few x 10 6 years, but with systematically younger ages (< 10 6 years) at the lowest 
masses and larger ages (~ 10 7 years) at the highest masses. The apparent mass dependent 
age inferred here is often seen in comparisons of cluster isochrones with pre-main sequence 
evolutionary models, and may be caused by incorrect assumptions in evolutionary models 
or the assumed temperature scale used for T Tauri stars. Finally, the temperatures and 
luminosities corresponding to the "stellar birthline" for spherical accretion (Stahler 1988) 
are shown, using the mass-radius relationship defined in Fletcher & Stahler (1994) as applied 
to the Siess et al. (2000) evolutionary models. 

The stellar luminosity estimates of the environmentally young stars are consistent, on 



2 The mean T Tauri luminosities does not include the stars HBC 358, 359, 360, 361, 362, 372, and 392, 
which have insufficient Lithium to be confidently considered T Tauri stars, and does not include the binaries 
HBC 354/355 or HBC 356/357, whose primaries lie below the main-sequence assuming a distance of 140 pc. 
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average, with the mean T Tauri star luminosities. At masses below ~ 1 M Q , the stellar 
birthline is systematically over-luminous by 0.3 — 0.5 dex compared to the mean T Tauri 
luminosities and the environmentally young star luminosities, though statistically consistent 
with both of these distributions. At masses above ~ 1 M , however, the stellar birthline 
is systematically larger by 0.5 — 1.0 dex than the mean T Tauri luminosities, and only 
marginally consistent with the environmentally young luminosities. We note that the 10 6 
yr isochrone reasonably matches the mean luminosities of environmentally young stars over 
all masses. For direct comparison (and further use in Section 3.3 below), we list the 10 6 yr 
luminosity corresponding to the stellar temperature for all stars with temperature estimates 
in Table 4 (L]^ r ). Independent of selection criteria, environmentally young stars appear to 
be coeval with populations of Class II and III stars in Taurus. 



3.2.2. Stellar Masses and Angular Momentum 

The stellar luminosity and temperature estimates of the environmentally young stars 
allow direct comparison with pre-main sequence evolutionary models to estimate stellar 
masses. However, we do not favor computing masses via direct comparison given the large 
spread in stellar luminosities, likely caused by biases introduced from circumstellar material. 
The cool temperatures of these K and M type young stars suggest that most should reside 
on the convective and predominantly vertical (constant temperature) portion of their pre- 
main sequence evolution. Their temperature principally will determine their stellar mass. 
Therefore we use the temperature - mass relation defined by the Siess et al. (2000) 1 Myr 
isochrone to estimate masses for the stars in our sample with spectral type estimates. As 
noted above, the 10 6 yr isochrone reasonably represents the average apparent age of these 
stars. Masses for objects cooler than the lowest mass track (0.1 M ) are estimated to be 
0.07 M at spectral type M5.5 and 0.05 M at spectral type M6, based on both a modest 
extrapolation of the Siess et al. (2000) models and comparisons with other evolutionary 
models (e.g. Baraffe et al. 1998). 

The inferred stellar /substellar masses are listed in Table 3, and span nearly 2 orders of 
magnitude: 0.05 M to 3.5 M . The distribution of masses is shown in Figure 15 as a function 
of the two evolutionary diagnostics Tboi and a. For comparison, the additional sample of 
Class II stars is shown; their masses are determined from their assigned spectral types and 
the same temperature - mass relationship used above. The similar spectral type distribution 
of Class I and Class II stars (Section 3.1.2) translates into a similar mass distribution. 
Likewise, the slight trend towards hotter spectral types for HH stars translates into a slight 
difference in the stellar mass distributions; HH stars are slightly more massive. It should 
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be realized, however, that the comparison Class II population shown in Figure 15 is not 
representative of the true distribution of stellar (and substellar) masses in Taurus. The 
now well established populations of low mass stars and brown dwarfs in Taurus (Briceno 
et al. 1998, 2002; Martin et al. 2001; Luhman et al. 2003, e.g.) are in general too 
faint to have been detected in previous mid- and far- infrared surveys (e.g. IRAS), from 
which spectral indices and bolometric temperatures can be determined. Thus without these 
evolutionary diagnostic measurements, they can not be plotted in Figure 15. Nevertheless, 
the 3 Class I stars IRAS 04158+2805, IRAS 04248+2612, and IRAS 04489+3042 are of 
particular interest. With spectral types of M5.5 or cooler, they have substellar masses and 
are the first spectroscopically confirmed Class I brown dwarfs. 

Figure 16 shows the vsini values measured in Section 3.1.1 normalized by their break-up 
velocity, defined as v br = ^GM/R. Stellar mass and radius estimates are from the adopted 
Siess et al. (2000) 10 6 yr isochrone described above. The normalized velocities range from 
less than 5% to 20% of the break-up velocities. K-S comparisons of Class I stars with Class II 
stars, and HH stars with non-HH stars show no differences in distributions (K-S probabilities 
of 0.079 and 0.096, respectively). The difference in the vsini distributions of HH stars and 
non-HH stars (Section 3.1.1) disappears when these values are normalized by the break-up 
velocity, possibly because of the slight mass difference between HH and non-HH stars. There 
is no evidence for a change in the angular momentum between the evironmentally young and 
the T Tauri samples. 

3.3. Inferred Circumstellar Properties 

In this section mass accretion rates are estimated from the observed continuum excesses 
and mass outflow rates are estimated from the strengths of forbidden-line emission, under 
the assumption of a bipolar flow. We outline the steps of these procedures in some detail to 
highlight the substantial uncertainties in current estimates. The ratio of these quantities are 
used to constrain the origin and energetics of these apparently physically related processes. 

3.3.1. Mass Accretion Rates 

The continuum excesses observed in T Tauri star spectra are attributed to the high- 
temperature regions generated as accreting material shock decelerates at the stellar surface. 
Measurements of this excess can therefore be used to estimate the mass accretion rate under 
the assumption that the liberated energy is gravitational potential energy. The steps involved 
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in estimating mass accretion rates are illustrated in the 5 panels of Figure 17. Figure 17 
contains only a subset of our data; it is restricted to a narrower mass range for direct 
comparison with the stellar and accretion properties of the 17 Taurus T Tauri stars studies 
by Gullbring et al. (1998). The mean and range from this study are indicated in each panel. 

The top panel of Figure 17 shows the distribution of continuum excesses at 6500 A 
( r 650o) as a function of stellar mass for stars observed here and the additional sample of 
Class II stars 3 (Section 3). The distribution of continuum excesses is indistinguishable from 
that measured by Gullbring et al. (1998) at bluer wavelengths (3200 - 5300 A). As noted 
by (Basri & Batalha 1990) and further supported by the results in Section 3.1.2, the veiling 
is relatively constant redward of ~ 5000 A. 

Since the continuum excess measurements are determined relative to the photosphere 
(r = F excess I Fphotosphere) , the photospheric flux needs to be known in order to determine the 
continuum excess flux. Since current estimates of the photospheric flux are subject to very 
large uncertainties (Section 3.2.1), we assign stellar flux values based on the predictions of 
the Siess et al. (2000) evolutionary models for 1 Myr aged stars as observed through the R c 
passband. An R c zero point of 2.32 x 10~ 9 erg s _1 cm -1 A -1 (Hayes 1985) is assumed. The 
second panel of Figure 17 shows the distribution of stellar luminosities corresponding to the 
adopted photospheric flux levels. We show stellar luminosities, rather than i? c -band fluxes, 
for direct comparison with the Gullbring et al. (1998) values. Our luminosities are slightly 
larger than those determined by Gullbring et al. (1998), but by less than a factor of 2 in the 
mean. The third panel in Figure 17 shows the distribution of excess fluxes measured over 
the 6000 - 6500 A wavelength interval. Surprisingly, the distribution is indistinguishable 
from that measured over the 3200 - 5300 A interval by Gullbring et al. (1998); standard 
hot spot models predict the excess flux should be only one-tenth of that seen in the bluer 
passband (Hartigan & Kenyon 2003). Up to a factor of 2 of this discrepancy may be caused 
by larger stellar luminosities, but the observed red excesses are still larger by a factor of ~ 5 
compared with standard hot spot model predictions. 

In order to convert the excess measurements to the total accretion luminosity, a bolo- 
metric correction is needed. Based on the calculations of Hartigan et al. (1991) who model 



3 The continuum excess measurements for the additional Class II stars are from Basri & Batalha (1990); 
Hartigan et al. (1995) and White & Basri (2003). Although the values measured by Basri & Batalha 
(1990) and White & Basri (2003) are determined near 6500 A, as are those newly presented here, the values 
measured by Hartigan et al. (1995) are determined at a slightly shorter wavelength, 5700 A. However, direct 
comparison of the veiling values for stars observed by both Hartigan et al. (1995) and Basri & Batalha 
(1990) shows no difference (16 star overlap; median r 6 5oo - ^5700 = +0.03). We use the rs7oo values as r65oo 
values. 
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the emission from the accretion shock as a slab of pure hydrogen gas of constant temperature 
and density in LTE, Hartigan & Kenyon (2003) estimate that the bolometric correction for 
an excess measured over 6000 - 6500 A is ~ 35. However, this bolometric correction, which 
is 10 times larger than the bolometric correction predicted for the 3200 - 5300 A passband, 
will yield accretion luminosities that are 10 times larger than that measured by Gullbring 
et al. (1998) since the two passbands have similar flux excesses. Consequently, we suspect 
that the mass accretion rates determined by Hartigan & Kenyon (2003) are most likely 
over-estimates, as their best matched model would predict very large blue wavelength ex- 
cesses, in contrast to what is observed. Similarly, we suspect that inferred mass accretion 
rates determined by Gullbring et al. (1998) are most likely under-estimates, as their best 
matched model predict negligible red excesses, in contrast to what is observed. Without 
a better model for the total continuum excess flux, we adopt bolometric correction of 11, 
corresponding to the logarithm average of the two extremes (3.5 and 35). The resulting 
accretion luminosities are shown in the fourth panel of Figure 17 and are roughly a factor of 
3 larger than those estimated by Gullbring et al. (1998). 

Finally, the accretion luminosities are converted to mass accretion rates (Figure 17; 
bottom panel) by assuming that the accretion shock luminosity equals the gravitational 
energy released per second as material free-falls along magnetic field lines from an inner disk 
radius of 3R sta r (e.g. Gullbring et al. 1998; Hartigan & Kenyon 2003). This assumption will 
under-estimate the accretion rate if a non-negligible fraction of the gravitational energy is 
used to power a wind or jet (the accretion rate will increase by a factor of 2 in the case where 
1/2 the energy is used in mass loss processes). Comparison of the mechanical luminosity 
of the outflow to the accretion luminosity suggests this is unlikely (Section 3.3.3). The 
stellar radii used in these calculations are determined from the temperature - radii relation 
of the Siess et al. (2000) 10 6 yr isochrone. Our geometric assumptions differ from those 
of Gullbring et al. (1998) only in the choice of inner disk radius. This small difference, in 
combination with the difference in accretion luminosities, leads to mass accretion rates that 
are, on average, a factor of 4 larger than those of Gullbring et al. (1998). The fifth panel of 
Figure 17 shows the distribution of inferred mass accretion rates. 

The mass accretion rates for stars observed here are listed in Table 4 and span over 
3 orders of magnitude, from 4 x 10~ 10 M /yr to 7 x 10~ 7 M /yr. The distribution of 
mass accretion rates determined for environmentally young stars and for the additional 
sample of Class II stars are shown in Figure 18 as a function of the evolutionary diagnostics 
T bo i and a. The median mass accretion rate of Class I stars is slightly less than that of 
Class II stars (median Class I log(M4 CC ) = -8.1; median Class II log(M4ce) = -7.3). This 
is primarily because of the strong mass dependence of the accretion rate (e.g. White & 
Basri 2003; Muzerolle et al. 2003); unlike the Class I sample, the Class II sample does 
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not extend into the substellar regime because of insufficient data for determining bolometric 
temperatures and spectral indices. Excluding the 3 Class 1 brown dwarfs, the distributions of 
mass accretion rates for Class I stars and Class II stars are indistinguishable (K-S probability 
of 0.60) and have more similar median values (median stellar Class I log(M^ cc ) = -7.1). The 
distributions of mass accretion rates for HH stars and non-HH stars are also similar (K-S 
probability of 0.11) though HH stars have mass accretion rates that are larger by a factor of 
2.5 in the median (median HH log(M j4cc ) = -7.0; median non-HH \og(M Acc ) = -7.4). Table 3 
summarizes the median statistics, excluding brown dwarfs. Combining all types, the typical 
(median) mass accretion rate for an optically veiled K7-M1 star is 4 x 10~ 8 M /yr, with a 
range of 2 orders of magnitude. 

The inconsistencies between the predictions of standard hot spot models and the large 
red excesses observed merit some discussion. Previous work has demonstrated that the 
variations in the veiling at blue wavelength are simultaneously correlated with variations at 
red wavelengths (see Figure 4 in Basri & Batalha 1990), suggesting a physically related 
origin. Emission from the inner disk is unlikely to contribute much optical excess flux if the 
temperature is limited to the dust destruction temperature of silicate-type grains (~ 1500 
K), as some observations suggest (Muzerolle et al. 2003b). The interpretation that we 
favor is that the this emission originates in a cooler component of the accretion generated 
shock. If this cooler component dominates the excess at red wavelengths, then the observed 
constant flux over the wavelength range 6500 A- 8400 A (Figure 10) requires a surprisingly 
cool temperature of only ~ 4000 K. Thus for a K7 star of similar temperature, the excess 
emission should be constant with wavelength, as observed (Figure 10; Basri & Batalha 1990). 
This cooler component must be featureless, suggesting that it is heated from above. The 
cooler excess component must also cover a large faction of the star in order to be observed. 
For example, in the case where the star and the excess have the same effective temperature 
(e.g. for a K7 star), the veiling, or the ratio of excess flux to photospheric flux, will be 
equal the ratio of the area of these emitting regions (r = F ex /F phot = A ex /A phot ). Thus for 
r = 1.0, the cool featureless emission will cover 1/2 of the stellar surface. A more detailed 
investigation of the accretion generated excess will require broader wavelength coverage that 
what is presented here (e.g. Stempels & Piskunov 2003) and is consequently beyond the 
scope of this study. 



3.3.2. Mass Loss Rates 



Optically thin forbidden-line emission can be used to estimate the mass loss rate from 
young accreting stars. Forbidden emission lines such as [OI] 6300 A and [SII] 6731 A in 
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the spectra of accreting T Tauri stars nearly always show emission at the systemic velocity 
and often show a separate blue-shifted peak or blue-shifted asymmetry (corresponding to 
velocities of a few x 10 to a few x 100 km/s; see Appenzeller et al. 1984; Mundt et al. 
1987; Edwards et al. 1987; Cabrit et al. 1990; Hamann 1994; Hartigan et al. 1995). 
The blue shifted emission is attributed to a high velocity wind or jet on the near side of 
the star; the presence of a circumstellar disk is believed to prevent direct observation of 
the redshifted component of the presumed bipolar flow. High spatial resolution imaging 
studies of young stars with strong forbidden-line emission have shown that the high and 
low velocity components likely originate in separate regions. The low velocity emission is 
usually spatially coincident with the stellar position, while the high velocity emission appears 
spatially extended, often in the form of a well collimated jet (e.g. Mundt et al. 1987; Lavalley 
et al. 1997; Bacciotti et al. 2000). The presence of a high velocity component of course 
requires the flow to be directed nearly along our line of sight; several systems with edge- 
on disk orientation power optical jets but have no high velocity forbidden line emission, as 
expected (Section 3.1.4). Interpretation of the low velocity forbidden line-emission is still 
uncertain (see e.g. Hartigan et al. 1995); it too may originate from outflowing material at 
slower velocities, possibly a disk wind (Kwan & Tademaru 1995). 

We estimate mass outflow rates using the prescription given by Hartigan et al. (1995) for 
[SII] 6731 A emission. In our spectra, this forbidden emission-line is typically the strongest 
and the least contaminated by telluric or stellar features. Since the [SII] emission is op- 
tically thin, the observed [SII] luminosity can be used to estimate the total mass of the 
emission region. Hartigan et al. (1995) predict the mass of the emitting region to be 
1.43 x 10~ 3 (L673i/-L©) M Q , where L^^i is the [SII] 6731 A luminosity in solar units. This 
mass estimate assumes the electron density is in the high density limit (N e > 2 x 10 4 cm~ 3 ; 
consistent with [SII] line ratios), that all the Sulfur atoms are singly ionized, and a cosmic 
sulfur abundance. With a mass M of the emitting region, the mass loss rate can then be 
estimated from the flow speed V and the length scale L over which the emission is observed 
(Mont = MV/L). 

The EW[SII]s are converted to fluxes by assuming that the underlying continuum flux 
level is that of the star in the R c passband, as predicted by the Siess et al. (2000) 1 
Myr isochrone, multiplied by (l+r 650 o) to account for veiling. This is the same continuum 
assumption used to estimate the mass accretion rates. The majority of stars observed do 
not have well-separated high- and low-velocity forbidden-line emission, which would allow 
us to directly measure the velocity of the flow and the fraction of high velocity flux relative 
to the total line flux. Therefore for all sources we assume a flow speed of 150 km/s and 
a high velocity flux equal to 40% of the total [SII] 6731 A line flux. This flow speed is 
the same as that used by Hartigan et al. (1995) and is consistent with the de-projected 
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velocities of propagating HH knots (e.g. Krist et al. 1999; Bacciotti et al. 1999). The 
adopted percentage of high velocity flux is the mean percentage of high velocity emission 
measured by Hartigan et al. (1995, compare their Tables 3 and 4). The scale length is 
assumed to be the slit- width of our observations (1'.'15), projected to the distance of Taurus 
(140 pc). We note that the slit-width is of similar size to the seeing disk. Thus even if the 
jet was projected along the length of the slit, the characteristic length scale would still be 
appropriate since the width of the aperture extraction is defined by the seeing-limited width 
of the continuum, and is relatively unaffected by spatially extended emission lines. 

Following this simple prescription, mass outflow rates are determined for all stars with 
[SII] 6731 A emission, and upper limits to the mass outflow rate are determined in cases 
where only detection limits are available. We note that for edge-on systems, the apparent 
bias towards larger forbidden emission-line EWs will correspondingly bias the inferred mass 
outflow rates. The determined mass outflow rates are listed in Table 4 and range from 
2 x 10~ 10 M Q /yr to 3 x 10 -6 M Q /yr. Ignoring edge-on systems, the largest mass outflow 
rate is then 6 x 10 -7 M /yr (DG Tau). The calculated mass outflow rates are shown in 
Figure 19 as a function of and a. Likely edge-on systems are indicated. [SII] 6731 A 
measurements for the additional Class II star sample are taken from Hartigan et al. (1995); 
the outflow rates are determined similarly. We note that the inferred measured mass outflow 
rates are slightly larger than, by a factor of 3, the values determined from high spatial 
resolution spectroscopy (e.g. Bacciotti et al. 2002; Woitas et al. 2002). The mass outflow 
rates are also roughly a factor of 2 larger than those measured by Hartigan et al. (1995) 
using the high velocity component of [OI] 6300 A. The latter difference is primarily due 
to the slightly higher stellar luminosities of our adopted isochrone. Both samples show a 
very large dispersion in the mass outflow rate (~ 3 orders of magnitude for K7-M1 spectral 
types). Thus, even though we use a weaker forbidden emission line to measure outflow rates, 
the similar dispersions suggest that the large range of values are dominated by stochastic 
differences between systems. 

The median mass outflow rates of Class I stars and Class II stars are similar (median 
Class I log(M 0ut ) = -8.0; median Class II log(M 0ut ) = -8.2). Although the [SII] 6731 A 
emission lines, from which the mass outflow rates are calculated, are systematically larger 
for Class I stars than Class II stars, this dependence is less significant than the stellar 
mass dependence of the mass outflow rate. If the 3 Class I brown dwarfs are excluded, the 
median mass outflow rate of Class I stars is 31 times greater than that of Class II stars 
(median stellar Class I log(M OMt ) = -7.0). The distributions are nevertheless statistically 
indistinguishable (K-S probability of 0.29; edge-on disk systems have been excluded). The 
outflow rates of HH stars are, in general, greater than those of non-HH stars (median HH 
star log(Mout) = -7.4; median non-HH star log(Mout) = -8.7). These distributions show 
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a more significant difference (K-S probability of 0.00049). Table 3 summarizes the median 
statistics, excluding brown dwarfs. Combining all types except edge-on systems, the typical 
(median) mass outflow rate for an optically veiled K7-M1 star is 2 x 10~ 9 M Q /yr. 

3.3.3. Comparison of Mass Accretion and Mass Loss Rates 

The measured mass outflow rates are correlated with the measured mass accretion rates 
over the 3 orders of magnitude in Macc inferred here. This correlation has been demonstrated 
previously (e.g. Hartigan et al. 1995) and suggests a physically related origin, though it in 
part stems from the mass dependence of the accretion and outflow rates. This correlation 
is not uniquely one-to-one, however. There are 14 optically veiled stars which show no [SII] 
emission, and 3 stars with optical veiling upper limits which show strong [SII] emission. 
Hartigan & Kenyon (2003) note similar disparities using the stronger [OI] 6300 A emission 
line. 

Although both the inferred mass accretion and outflow rates depend directly on the 
photospheric luminosities, assumed in our analysis to correspond to 1 Myr aged stars, the 
ratio of these rates is independent of the photospheric luminosity. These mass flow ratios are 
listed following the mass accretion and mass outflow rates in Table 5. In cases where [SII] 
6731 A emission is measured, but only veiling upper limits are available, the ratio is a lower 
limit. In cases where the star is optically veiled, but no [SII] 6731 A emission is detected, the 
ratios are upper limits. These values are plotted versus the evolutionary diagnostics Tboi and 
a in Figure 20. The ratio Mo ut j Ma cc for 5 of the 6 edge-on disk systems are ~ 10. Since the 
majority of disk material is expected to accrete onto the star, as opposed to being ejected 
in an outflow, the large ratios for these systems support the hypothesis that the forbidden 
emission-lines are biased towards artificially large values. Excluding edge-on disk systems, 
the measured ratios range from 3.2 x 10~ 3 to 1.7. The median log ratio (including limits) is 
-1.2, corresponding to a ratio of 0.05. This value is similar to, but larger by a factor of ~ 5, 
the average value measured by Hartigan et al. (1995). 

Although the distributions of mass flow ratios for Class I stars and Class II stars span 
similar ranges, the mass flow ratios of Class I stars are larger, in the median, than those 
of Class II stars (median Class I \og(M 0u t/M Acc ) = -0.04; median Class II \og(M 0ut / M Acc ) 
= -1.2). The mass flow rates of HH stars are similarly larger than the values for non-HH 
stars (median HH star \og(M 0u t/ Ma cc ) = -0.13; median non-HH star log (M 0u t/ Ma cc ) = 
-1.4). These differences are correlated however, since many Class I stars are also HH stars. 
Excluding Class I stars, the median HH star \og(Mo u t/ 'Ma C c) decreases to -0.81. If many of 
the Class I stars observed here have edge-on orientations (Section 4.5), the effect of this on the 
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observed forbidden-line EWs (Section 3.1.4) could explain their higher ratios of M 0u t/ Ma cc - 

The ratio of M Qut / M Acc appears to be independent of the mass accretion rate. Over 
the M Acc range 10~ 8 — 10~ 6 M /yr, the mass flow ratios are uniformly distributed. We 
note that the large range in Mo u t / Ma cc is not because of observational error, but a con- 
sequence of significant dispersion in the relative strengths of the forbidden emission-lines 
and the continuum excesses. The two observations of IRAS 04303+2240, which measure 
log (Mout/ M Acc ) at -2.0 and -0.5, suggest this ratio changes for a given star-disk system 
(thereby excluding random orientation effects). These changes may arise from either actual 
changes in the mass accretion and/or mass outflow rates, or from changes in the properties 
of the emission regions (e.g. temperature and density). The highly variable mass accretion 
rate will cause an additional source of variation for any given epoch. The previously ejected, 
spatially extended forbidden emission should be more correlated with a time averaged mass 
accretion rate, as opposed to the instantaneous measured mass accretion rates we present 
here. 

Theoretical estimates for the ratio of mass outflow to mass accretion rate range from 
~ 1/3 (e.g. Shu et al. 1994) for a collimated jet powered by a magnet ospheric accretion flow 
to ~ 10 -4 for a slow disk wind from a magnetically threaded accretion disk (Wardle & Konigl 
1993, e.g.). Although the high velocity components and the current ranges of log(Mo u t/ Ma cc ) 
appear to exclude low velocity disk winds as a dominant mass loss mechanism, the large 
range in mass flow ratios currently offer little constraint on proposed models for powering 
the observed jets. 

Finally, the ratio Moutj Ma cc offers a useful measure of the mechanical luminosity of 
the outflow (l/2M OMt V 2 ) relative to the energy generated by accretion (L Ac c)- The right 
hand ordinate of Figure 20 shows this ratio, assuming accretion onto a 0.5 M Q star of radius 
2.0 R Q , and a mass loss velocity (V) of 150 km/s. The median luminosity ratio is 0.076. 
Thus the mechanical luminosity of the outflow is approximately 8% of the observed accretion 
luminosity. 

4. Discussion 

The insights gained from our spectroscopic observations address various long-standing 
assumptions regarding the earliest optically visible stages of star formation. In order to place 
our new results in context, we present a brief review of the standard paradigm of low mass 
star formation and the taxonomy of young stars. 

The formation of a low mass star (M < 3 M Q ) is thought to begin with the dynamical 
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inside-out collapse of a dense molecular cloud core (e.g. Shu et al. 1987). A 3 stage 
description of evolution from this point towards the main sequence was proposed to explain 
3 relatively distinct types of young stars as distinguished by their optical through millimeter 
SEDs (e.g. Lada 1987; Adams et al. 1987; Wilking et al. 1989). Young stars with SEDs 
that peak at far-infrared wavelengths (30 — 100/im) are presumed to be in the initial main 
accretion phase. These stars are often called Class I stars or protostars. The majority of 
their luminosity escapes at far infrared wavelengths as reprocessed radiation in the form of 
thermal dust emission; extinction to the star is often high. Young stars with SEDs that 
peak at near infrared wavelengths are Class II stars or classical T Tauri stars. Their SEDs 
typically resemble a revealed stellar photosphere with significant infrared and millimeter 
excess emission attributed to thermal emission from a circumstellar disk. Young stars with 
SEDs that peak in the optical, with little or no evidence for an infrared excess are called 
Class III stars or weak-lined T Tauri stars. These stars are presumed to have dispersed 
the majority of their circumstellar disk material. Subsequent to this 3 Class classification 
scenario, Andre et al. (1993) proposed a new class of less evolved objects, called Class 
stars. These are stars that are so deeply embedded that they are only observable at far- 
infrared and millimeter wavelengths. Distinction between Class and Class I stars is often 
observationally ambiguous, however (e.g. Young et al. 2003); both are believed to be in 
the main accretion phase, though the envelope is presumed to be more mass than the star 
during the Class stage. 

The association of strong outflow/jet signatures, such as HH flows, with many protostars 
suggest that this process plays an important role in clearing envelope material. Within the 
Class O/I/II/III classification scheme, HH flows first appear associated with Class stars 
(e.g. IRAS 04368+2557 Eiroa et al. 1994), are most common among Class I stars, are 
associated with the about 10% of Class II stars (Gomez et al. 1997), and are not associated 
with any Class III stars. HH stars are therefore expected to be, on average, younger than 
most Class II or Class III T Tauri stars. 

Since this classification scheme describing the evolution of circumstellar material around 
solar-type stars was first proposed, a wealth of direct observational evidence has been ob- 
tained in support of it. However, a correlation with stellar age has never been established. In 
this section we use the stellar properties (mass, rotation, age), and circumstellar properties 
(accretion rate, outflow rate, disk orientation, disk mass) of the 11 Class I stars and 43 Class 
II stars studied here to address issues related to the evolution from the protostellar stage to 
the T Tauri stage. We use these same properties to search for differences between HH stars 
and non- stars. 
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4.1. The Masses of Class I Stars 

If the inside-out collapse of a singular isothermal sphere proceeds under the control of 
slow ambipolar diffusion processes (e.g. Shu 1977; Shu et al. 1987), the resulting infall rates 
are set by the local isothermal sound speed and therefore should remain roughly constant in 
time throughout the cloud. One implication is that more massive stars require more time to 
form than less massive stars. Thus, assuming coeval star formation, these conditions would 
yield Class I stars that are more massive, on average, than Class II stars. However, recent 
observational, theoretical and numerical work suggests that supersonic turbulent flows rather 
than static magnetic fields control star formation (see the review by Mac Low & Klessen 
2004), affecting both the infall rate and the formation time scale. As an example, numerical 
models of gravoturbulent fragmentation suggest that the initial infall rates may be time 
dependent (initially large), scale with the final stellar mass, and yield a main accretion 
timescale that is nearly independent of the final stellar mass (Schmeja & Klessen 2004). 

Our results are more consistent with the proposal that supersonic turbulent flows control 
star formation, rather than static magnetic fields. As shown in Figure 15, the distributions 
of stellar mass for Class I and Class II stars are statistically indistinguishable. Although 
the sample of Class I stars with mass estimates determined here is small, there are 3 brown 
dwarfs, 6 sub-solar mass stars (0.2 - 0.9 M ), and 2 super-solar mass stars. Previous work has 
determined the masses of 2 Class I stars based on disk kinematics under the assumption of 
Keplerian rotation (0.2-0.4 M Q for IRAS 04381+2540 and 0.35-0.7 M Q for IRAS 04365+2535; 
Brown & Chandler 1999). This distribution of Class I star masses agrees very well with 
the initial mass function measured from T Tauri stars in Taurus (Briceno et al. 2002). The 
implication is that the known population of optically revealed T Tauri stars is likely to be a 
good representation of the final initial mass function produced by the cloud. The current and 
presumably subsequent generations of Class I stars are unlikely to skew the mass function. 

The cool spectral types of 3 Class I stars suggest that they are Class I brown dwarfs 
(IRAS 04248+2612, IRAS 04158+2805, and IRAS 04489+3042). With the exception of the 
most massive of these 3, IRAS 04248+2612, the 1.3 millimeter continuum fluxes suggest 
circumstellar disk+envelope masses of only a few hundreths of a solar mass (Motte & Andre 
2001). Thus, even if the majority of circumstellar material were to accrete onto these brown 
dwarfs, both IRAS 04158+2805 and IRAS 04489+3042 will remain substellar. If all of 
the 0.25 M Q envelope accretes onto IRAS 04248+2612, it will be a low mass star. IRAS 
04248+2612 is also of interest as it is the lowest mass object known to power a molecular 
outflow (Moriarty-Schieven et al. 1992) and a HH object (HH 31; Gomez et al. 1997). We 
suggest that the lowest mass Class I stars/brown dwarfs are low mass because they formed 
from lower mass cloud cores, and not because they are at an earlier evolutionary stage and 
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have not yet accreted the majority of their final mass. The implication is that the formation 
process of substellar mass objects, down to ~ 0.05 M Q , is simply a scaled-down version of 
that for solar mass stars. Brown dwarfs likely form via the dynamical collapse of a cloud 
core, experience an embedded Class I phase of evolution, and are capable of powering both 
molecular and HH flows. The presence of circumstellar material and accretion argues against 
low mass star/brown dwarf formation scenarios involving early ejection from a competitively 
accreting cluster (e.g. Reipurth & Clarke 2001). 

Finally we note that without stellar temperature estimates, previous attempts to es- 
timate the stellar masses of protostars have proceeded by assuming accretion dominated 
luminosities; the stellar mass can be estimated by calculating the gravitational potential 
well required to liberate the observed luminosity. The low luminosities of Class I stars in 
Taurus consequently lead to very low substellar masses (e.g. Haro 6-33 = IRAS 04385+2550; 
Young et al. 2003, estimate M= 0.01 M ). These masses are generally inconsistent with 
the values determined more directly here (e.g. Table 4; for Haro 6-33, M= 0.46 M Q ). Mass 
estimates assuming accretion dominated luminosities appear to be inappropriate for many 
Taurus protostars. 



4.2. The Ages of Class I Stars 

Unlike the self-embedded Class I stars, the optically revealed nature of Class II and 
Class III stars have allowed more direct measurements of their stellar temperatures and 
luminosities for comparison with the predictions of pre-main sequence evolutionary models. 
These comparisons imply T Tauri ages of a fewxlO 6 yr, with no statistically significant 
difference between the HR diagram locations of Class II and Class III stars (Kenyon & 
Hartmann 1995; White & Ghez 2001; Briceno et al. 2002). Class II and III stars also 
have similar lithium abundances (e.g. Strom et al. 1989; Martin et al. 1994) supporting 
similar ages. Given that Class I stars in Taurus are approximately one-tenth as numerous as 
T Tauri stars (both Class II and Class III stars combined), it has been postulated that the 
Class I phase lasts for roughly 10% of the T Tauri phase, under the assumption of a constant 
star formation rate (e.g. Myers et al. 1987; Kenyon et al. 1990). The implied statistical 
age of Class I stars is then a fewxlO 5 yr. 

With stellar properties now determined for 42% of the Class I stars in Taurus, their 
stellar ages can be assessed more directly. As illustrated in Figure 14, environmentally 
young stars have ages that scatter uniformly about those of Class II and Class III stars, 
independent of the selection criteria used to identify them (Tboi, a, HH properties). The 
mean stellar luminosities of Class I stars, specifically, are also more consistent with those of 
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optically revealed T Tauri stars than with the stellar birthline, though the difference is only 
signficant above 1.0 M . Class I stars appear to be coeval with Class II and Class III stars. 

4.3. The Early Evolution of Angular Momentum 

The accretion of high angular momentum material from a rotating, collapsing protostel- 
lar envelope during the earliest stages of star formation is expected to produce very rapidly 
rotating stars (e.g. Durisen et al. 1989), with rotational velocities comparable to the break- 
up velocity (vb r = a/ GM/R ~ 300 km/s). Thus, it came as a surprise to discover that most 
T Tauri stars, at least in the Taurus star forming region, are slowly rotating. Nearly all of 
these stars have rotational velocities less than one-tenth tv (e.g. Hartmann et al. 1986; 
Bouvier 1990; Bouvier et al. 1995). Proposed explanations for the slow rotation rates 
usually involve "magnetic braking" , a mechanism in which the star is magnetically coupled 
to the accretion disk and transfers angular momentum to the more slowly rotating parts of 
the outer disk (Blandford & Payne 1982; Pudritz & Norman 1986), the inner disk (Konigl 
1991; Shu et al. 1994) or possibly a stellar wind (Kwan & Tademaru 1988). After the disk 
dissipates, substantial angular momentum loss must still occur, likely through a magnetically 
coupled stellar wind, in order to explain the slow rotation of somewhat older stars such as 
those in the Pleiades and Hyades (Stauffer et al. 1997; Krishnamurthi et al. 1998). 

Thusfar, the rotational properties of stars during the main accretion phase have been 
difficult to ascertain. High dispersion near-infrared spectra of stars in the p Ophiuchus dark 
cloud have found that at least some environmentally young stars rotate moderately rapidly 
(t>sim = 30-50 km/s; Greene & Lada 1997, 2000), with tentative evidence that stars with 
large spectral indices rotate the fastest (Doppmann et al. 2003), though this is based on only 
3 stars with spectral indices > 0.0 (i.e. Class I-like). Here we have demonstrated that Class 
I stars in Taurus are, as a group, slowly rotating, with a i>sim distribution indistinguishable 
from that of the more environmentally evolved Class II stars (Figure 9). None of the observed 
Class I stars in Taurus rotate with vsini >30 km/s. Further, since the Class I luminosities 
and hence radii are similar to those of the Class II stars, there is no evidence for an change 
in angular momentum from the Class I to the Class II stages (e.g. Figure 16). The youngest 
optically visible stars have already been slowed to well below (<5-15%) break-up velocity. 
Stars in Taurus are known to be systematically slower rotators than stars in Orion (Clarke & 
Bouvier 2000; White & Basri 2003). Similar differences in the cloud properties or possibly 
stellar ages (see below) may explain apparent differences in the rotation rates of Taurus Class 
I stars and p Ophiuchus Class I stars, once confirmed with a larger p Ophiuchus sample. 

It has been proposed that T Tauri stars in Taurus rotate more slowly than stars in 
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other star forming regions (e.g. p Ophiuchus, Orion) because Taurus stars are older (by a 
factor of ~ 3), having ages larger than the expected disk braking timescale (Hartmann et al. 
2002). In this scenario, the slow rotational velocities of Class I stars therefore corroborates 
the suggestion in Section 4.2 that these stars are as old as T Tauri stars. We speculate 
that if these Class I stars indeed represent the longest-lived disk population of Class II stars 
(Section 4.6), this already slowly rotating sample could produce the "slow rotators" observed 
at intermediate pre- main sequence ages (see Bouvier 1997). 

4.4. The "Luminosity Problem" of Class I Stars, Revisited 

The mass infall rate during the main accretion phase can be approximated by dividing 
the local Jeans mass by the free-fall timescale (Mj n f a u ~ Mj/rjf = 5Ac 3 JG; Schmeja & 
Klessen 2004). For an isothermal sound speed, c s , of 0.2 km/s characteristic of Taurus, 
this corresponds to M In f aU = 1 x 10~ 5 M Q /yr. The effects of magnetic support can produce 
smaller infall rates (2 x 10~ 6 M /yr; Shu 1977) while the external compression in turbulent 
flows can produce larger infall rates (up to ~ 10~ 4 M /yr, at least initially; Schmeja & 
Klessen 2004). Thus if Class I stars are in the main accretion phase, they should have mass 
infall rates consistent with these predictions. In support of this, density profiles correspond- 
ing to infall rates of a fewxlO -6 M /yr in envelope-only models successfully fit the SEDs 
and scattered light images of many Class I stars (Kenyon et al. 1993a,b; Whitney et al. 
1997). 

As first pointed out by Kenyon et al. (1990), if the infalling envelope material is 
channeled onto the star via steady-state disk accretion at this rate, the liberated accretion 
luminosity would be roughly 10 times that emitted from the photosphere. With accretion 
dominated luminosities, Class I stars would have total luminosities roughly 10 times those 
of similar mass T Tauri stars. Direct tests of this prediction have been inhibited, however, 
by the lack of measured stellar and accretion luminosities of Class I stars. Nevertheless, 
previous studies have tried to assess the relative contribution of the accretion luminosity, 
under the assumption of similar stellar luminosities, by comparing the bolometric luminosi- 
ties of Class I stars and Class II/III T Tauri stars. Bolometric luminosities are calculated 
by integrating the entire energy distribution of a star. These comparisons show no statis- 
tically significant difference between these Classes. The discrepancy between the observed 
bolometric luminosities and the predicted accretion dominated luminosities is often referred 
to as the "Luminosity Problem" for Class I stars (Kenyon et al. 1990, 1994). 

The spectroscopic measurements of 11 Class I stars observed here allow a more direct 
measure of the relative stellar and accretion contributions to the bolometric luminosity. We 
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first confirm that the stellar, accretion, and disk emission can account for the observed 
bolometric luminosity. The total luminosity, L tot , is calculated as 

Ltot L 'star L accs fo ock L acc d<i sk ^rep^diski (1 ) 

where L star is the stellar luminosity calculated in Section 3.2.1 (Table 4), L acc ^h ock is the 
accretion luminosity calculated in Section 3.3.1 (L acC;Shock = GMM star [l/ R star — 1/R in }), 
L acc ,disk is the viscously generated luminosity of the disk, and L rep ^i S k is the reprocessed 
luminosity from the disk. The 3rd and 4rth terms are calculated assuming a geometrically 
thin, optically thick circumstellar disk with an inner radius (Ri n ) of 3 R s t ar and infinite outer 
radius. The luminosity viscously generated in the disk equals the change in potential energy 
of material moving from infinity to the inner disk edge, and thus L acc ^ isk = GMM star (1/R in — 
l/i?oo), which simplifies to L acc4isk = 0.5 x L acCtShock . Hence, L acc4isk can be determined 
from the inferred L acCyS h oc k values. The luminosity reprocessed by the disk originates from 
two sources, the central star and the bright accretion shock, and thus L rev ^ sk = 0.25(L star + 
L a cc,shock) Rstar / Rin (Adams & Shu 1986). The total luminosity can then be simplified as 

Ltot — l-08L star + 1.58L accs /j OC fc. (2) 

The total luminosities calculated following this prescription agree well with the bolo- 
metric luminosities for both Class I stars and Class II stars (median Class I log(L tot / 'L^) 
= —0.08, a = 0.58; median Class II log(L tot /L bo i) = 0.19, er = 0.43). The agreement suggests 
that the dominant luminosity sources are accounted for; the large scatter suggests that the 
absolute values of the luminosity sources are very uncertain. For the Class I stars, this agree- 
ment suggests there is not significant contribution to the bolometric luminosity generated 
from an infalling envelope accreting directly onto the disk at large radii. 

For both Class I and Class II stars, the contribution of disk accretion (L acCyShock + 
L aC c,disk) to the bolometric luminosity can range from a few percent to 50 percent. However, 
in the typical (median) case, only 25% of the bolometric luminosity is generated through 
disk accretion; the majority of the bolometric luminosity originates from the star (see also 
Muzerolle et al. 1998). We note that the bolometric luminosities of the 11 Class I stars 
studied here (median \og(Lboi/ L Q ) = — 0.36,(7 = 0.61) are very similar to the bolometric 
luminosities of the remaining 15 Class I stars in Taurus (median log(Lb o ;/L ) = — 0.17, a = 
0.56). Thus, if the unobserved Class I stars have a similar distribution of stellar luminosities, 
then they can not have disk accretion rates as high as predicted by envelope infall models 
either. These results strongly favor one proposed resolution to the luminosity problem - 
Class I stars do not have accretion dominated luminosities. However, this resolution raises 
2 additional questions: (1) are known Class I stars properly classified?, (2) are Class I stars 
in the main accretion phase? We address these 2 issues in the next sections. 
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4.5. The Effect of Orientation on Class Classification 

In the traditional classification scheme, Class I stars are true protostars - stellar embryos 
surrounded by an infalling envelope. However, the observational criteria traditionally used to 
identify Class I stars distinctly from Class II stars depend upon the disk/envelope orientation 
relative to line of sight. Radiative transfer models of still-forming stars predict that edge- 
on systems will have optical, near- and mid-infrared characteristics similar to less evolved, 
more embedded stars (Kenyon et al. 1993a,b; Yorke et al. 1993; Sonnhalter et al. 1995; 
Men'shchikov & Henning 1997; Whitney et al. 2003b). As examples, Class II stars viewed 
edge-on will, in many cases, resemble Class I stars, while Class I stars viewed at pole- 
on orientations will have optical/infrared properties more characteristic of Class II stars. 
The predictions for pole-on orientations are less clear, however, since the distribution and 
radiative properties of material from bipolar flows is not included in most models. We note 
that if Class II stars with disk inclinations of 90 ± 5° appear "edge-on", then 8.7% of all 
disk systems should appear edge-on, assuming random orientations (17.4% if % — 90 ± 10°). 
It is important to recognize how these orientation effects could have introduced selection 
biases into both the optically revealed Class I stars studied here, and the known population 
of Class I stars in Taurus. 

As discussed in Section 3.1.4, the emission-line profiles offer an indirect constraint on 
the star/disk orientation. The general trend of narrow emission line profiles at the systemic 
velocity suggests a bias towards edge-on orientations for the 11 Class I stars with measured 
photospheric features. Only 3 of these stars show high velocity (> 60 km/s) forbidden line 
emission, compared to more than half of Class II stars (e.g. Hartigan et al. 1995). The 
Ha 10%-widths of these Class I stars are systematically less than those of Class II stars 
(Section 3.1.4). Only 1 of the 11 Class I stars (HL Tau) shows evidence of a strong blue 
shifted absorption superimposed upon the Ha emission profile. In comparison, roughly half 
of Class II stars show this absorption feature (e.g. Alencar & Basri 2000). These combined 
results imply little or no high velocity material along our line of sight for most Class I 
stars, as expected for an edge-on orientation. A orientation bias could also explain the 
unusually large forbidden emission-line EWs of most Class I stars; the somewhat extended 
forbidden emission line region may be more directly observable in this case than the obscured 
stellar photosphere. Overall, the emission line profiles of most optically revealed Class I 
stars favor nearly edge-on orientations. One clear exception is HL Tau. Its high velocity 
forbidden emission lines, broad Ha profile with superimposed absorption favor a less edge- 
on orientation. This orientation may explain why this is one of the most optically bright 
environmentally young stars. 

The evidence for edge-on orientations does not necessarily imply that the observed 
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sample of Class I stars are in fact Class II stars. The presence of spatially extended envelope 
material, as determined from image morphology at infrared and millimeter wavelengths, 
offers a more direct constraint on the evolutionary Class. Based on criteria put forth by 
Motte & Andre (2001), only 58% (15/26) of the Class I stars listed in Table 1 are true 
protostars. The remaining 42% (11 stars) have envelope masses < 0.1 M and are spatially 
unresolved at 1.3-mm wavelengths (referred to as "unresolved Class I sources" in Motte & 
Andre 2001, these stars are marked as Class I' stars in Table 1). Motte & Andre (2001) 
suggest that these stars are more likely transitional Class I/II stars or highly reddened Class 
II stars. The complementary near-infrared morphology survey by Park & Kenyon (2002) 
supports the claim that these 11 stars are not bona fide Class I stars 4 . However, we note that 
the morphological criteria used in these studies do not account for the luminosity and mass 
of the central star. IRAS 04158+2805 and IRAS 04489+3042 may appear more evolved and 
point-like because they are lower luminosity Class I brown dwarfs with smaller disks and 
envelopes. 

We conclude that as many as one-half of the Class I sample listed in Table 1, and 
in particular 8 of 11 Class I stars studied here, could be misclassified Class II stars. The 
emission-line profiles and image morphology suggests that in many cases this misclassification 
may have been caused by a nearly edge-on orientation. We nevertheless continue to refer 
to all of these as Class I stars since the orientations and the mass effects are generally not 
known. 



4.6. Are Class I Stars in the Main Accretion Phase? 

Although the absolute values of the mass accretion rates remain considerably uncertain, 
typical Class I and Class II stars of K7-M1 spectral type have Macc ~ 4 x 10~ 8 M /yr 
(Section 3.3.1). There is no difference between the disk accretion rates of Class I stars 
and Class II stars. Muzerolle et al. (1998) found similar results based on Br7 luminosity 
measurements and assumed stellar properties. Although the new disk accretion rates we 
determine in Section 3.3.1 are somewhat larger than previous estimates (e.g. Gullbring et al. 
1998), they are still 2 orders of magnitude less than the envelope infall rates inferred from 
SED modeling (fewxl0~ 6 M Q /yr). 



4 Except for the 2 bona fide Class I stars with optically revealed companions (MHO 1/2 and GV Tau 
AB), which could bias T^oi estimates towards hotter temperatures, all bona fide Class I stars with T^oi 
estimates (11/13) have T^oi < 350 K. A T^oi < 350 K may be a better distinguishing value for Class I stars. 
A comparison of the spectral indices suggests no similar correction; the spectral indices of bona fide and 
non-bona fide Class I stars overlap signficantly. 
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We explore possible explanations for reconciling this discrepancy. One possibility sug- 
gested by Kenyon et al. (1990) is that the infalling envelope material is not transferred 
to the star via disk accretion in a steady-state fashion. The accreting envelope mass ac- 
cumulates in the circumstellar disk until it becomes gravitationally unstable (e.g. Larson 
1984) leading to a marked increase in the mass accretion rate. This process may be related 
to the FU Orionis outburst phenomenon, which is currently understood as rapid (~ 10~ 5 
M /yr), but short lived (~ 100 yr) accretion through a circumstellar disk onto a central 
star (Hartmann & Kenyon 1987). At this rate of accretion, Class I stars would only need 
to spend a 5-10% percent of their lifetime in the high accretion state to achieve typical T 
Tauri masses within 1 Myr. The most luminous Class I star in Taurus, L1551 IRS 5 (20.9 
L Q ; Kenyon & Hartmann 1995) has been proposed to be in an FU Ori-like outburst state, 
though a scenario in which it is an embedded star of a few solar masses satisfies its observed 
properties equally well. 

If the envelope material of Class I stars is accumulating in their circumstellar disks, 
they are expected to have more massive disks than Class II stars. One direct way to test 
this hypothesis is with 1.3-mm continuum observation. Since circumstellar dust is mostly 
optically thin at these wavelengths, the emission should trace the total mass in circumstellar 
dust. Figure 21 shows the distribution of 1.3-mm flux densities versus the evolutionary 
diagnostics and a for all stars in Taurus with these values, including the larger sample 
of Class I stars not spectroscopically studied here. 1.3-mm measurements of Class II stars 
are taken from Beckwith et al. (1990), Osterloh & Beckwith (1995), and Motte & Andre 
(2001), which all have beam sizes of ~ 11 — 12 arcseconds. The emission in this beam is 
likely to be a good tracer of circumstellar disk mass without significant contamination from 
any possible envelope mass that would dominate on larger spatial scales (Motte & Andre 
2001). Figure 21 also indicates (right hand ordinate) the corresponding dust mass for a dust 
opacity per unit mass column density of 0.01 cm 2 g _1 and a dust temperature of 20 K for all 
sources 5 . 

As Figure 21 clearly demonstrates, the distribution of circumstellar disk masses for Class 
I and Class II stars are indistinguishable. On average, Class I stars do not have more massive 
disks than Class II stars. One possible caveat is that if much of the disk mass of Class I stars 
is confined to be very close to the star (< 1 AU), as would be expected immediately prior 
to an FU Ori outburst, this mass would likely be optically thick and thus un-revealed in 
our 1.3-mm comparison. However, if the disks of Class I stars undergo FU Ori outbursts at 
semi-regular intervals (~ 10% of the time), the optically thin outer disks would still need to 



5 Although there is some evidence that both the dust opacity and temperature may change with evolu- 
tionary state, the expected changes are thought to be less than a factor of <~ 2 (Henning et al. 1995). 
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be more massive than those of Class II stars, on average, to sustain the larger time-averaged 
accretion rates. This is inconsistent with the observations and suggests Class I stars do not 
undergo FU Ori outbursts more often than Class II stars. 

An alternative possibility that we favor for reconciling the mass infall and disk accretion 
rates is that most of the Class I stars in Taurus are not in the main accretion phase. First, 
the current populations of Class I and Class II stars are likely biased because the available 
low spatial resolution mid-infrared measurements (e.g. IRAS) used in Class classification 
criteria (e.g. a, T& j) often include several young stars (Section 2.1). Some Class I and 
Class II stars may have been misclassified as a consequence. Additionally, as discussed in 
Section 4.5, 42% of stars classified as Class I stars do not appear to be bona fide proto- 
stars (ie. stars with spatially extended envelope structures). Many of these are likely to 
be Class II or transitional Class I/II stars that have been misclassified because of nearly 
edge-on disk orientations. Stars like IRAS 04016+2610 and IRAS 04302+2247, as examples, 
have morphologies and kinematics that are better described by a rotating disk-like structure 
(Hogerheijde & Sandell 2000; Boogert et al. 2002; Wolf 2003) than a collapsing envelope 
model (Kenyon et al. 1993a; Whitney et al. 1997). In the more general case, we suspect 
that the envelope infall rates have been over-estimated, even for bona fide Class I stars, by 
improperly accounting for the emission from a circumstellar disk. For an infalling envelope 

1 /2 

model, the ratio M In f a u/Mj sets the peak wavelength of the SED (Kenyon et al. 1993a); 
higher infall rates produce redder SEDs. However, the addition of an optically thick circum- 
stellar disk also shifts the peak of the SED toward redder wavelengths (cf. Kenyon et al. 
1993a; Wolf 2003). Thus, if most Class I stars have circumstellar disks as expected, then 
evelope-only models (ie. no circumstellar disk) will systematically over-estimate the mass 
infall rate in order to account for their redder SED. 

We emphasize that even if, as we propose, most Class I stars are not in the main accretion 
phase, they nevertheless must still have circumstellar envelopes. Disks alone are insufficient 
to explain the high extinctions, scattered light images, 1.3-mm morphologies, and mid- to far- 
infrared SEDs. Our proposal is simply that there is less mass in the envelope than has been 
presumed for many of the Class I stars in Taurus. Consequently, the mass infall rates have 
been over-estimated. Proper determination of the envelope mass and infall rate is especially 
challenging for the low luminosity Class I stars in Taurus, whose envelope emission can be 
confused with diffuse cloud emission (e.g. MHO 1/2, IRAS 04361+2547, IRAS 04381+2540; 
Motte & Andre 2001). Additionally, the large fraction of companion stars with ~ 10 3 AU 
separations (Duchene et al. 2004) can give the impression of spatially extended structures. A 
convincing case for a massive infalling envelope can only be established by spatially mapping 
molecular line profiles and accounting for the effects of outflows and rotations (Evans 1999). 
Of the stars listed in Table 1, only the Class star IRAS 04368+2557 (L1527) retains 
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a massive extended envelope and shows unambiguous evidence for infall (Gregersen et al. 
1997). 

If Class I stars are in fact as old as most T Tauri stars (~ 1 Myr; Section 4.1), the 
presence of envelope material in these systems requires a long envelope dispersal timescale 
(ie. comparable to T Tauri ages). A large dispersion in the disk dispersal timescale is thought 
to explain the apparently coeval populations of Class II and Class III stars. Extrapolating 
this idea to include high-latitude envelope material, many Class I stars may represent a T 
Tauri star sub-sample with the longest envelope/disk dispersal timescale. The reason for 
this large dispersion is unclear, but may involve variations in the initial properties of the 
cloud core, dynamical and radiative effects from the presence of close stellar companions, or 
the onset of planet formation. 

The combined properties of Class I stars in Taurus support a scenario in which most 
are not in the main accretion phase. For these systems, the stellar (or substellar) masses 
must have accumulated prior to the Class I phase, during the Class phase. The flat density 
profiles of pre-stellar cores (Ward-Thompson et al. 1994) suggest initial conditions that favor 
more rapid accretion during the Class phase (Henriksen et al. 1997). The stronger mass 
outflows of Class stars relative to Class I stars supports this interpretation (Bontemps et al. 
1996), if powered by mass accretion as suspected. A possible caveat to this suggestion is that 
if the Class phase is the main phase of accretion, the small number of known Class stars 
in Taurus suggests a very rapid formation timescale (~ 10 4 yr) and consequently, a very high 
mass accretion rate (~ 10~ 4 M /yr). These predictions are inconsistent with the relatively 
low luminosities of known Class stars (assuming they too are properly classified); Class 
stars would have their own luminosity problem. However, the IRAS satelite may not have 
been sensitive enough to detect most Class stars in Taurus (e.g. IRAM 04191+1522; Andre 
et al. 1999). Indeed, early Spitzer Space Telescope results show that some "starless cores" 
may not be not be starless (Young et al. 2004). The 44 starless cores in Taurus (Onishi et 
al. 2002) may contain a yet hidden population of Class stars. Thus, concerns about the 
implied statistical ages and infall rates of Class stars in Taurus may be premature. 



4.7. Herbig-Haro Energy Sources 

Many of the Class I and Class II stars studied here have been identified as the energy 
sources of HH objects (Section 2.1). Some studies suggest HH stars have higher mass accre- 
tion rates and larger circumstellar reservoirs than the average T Tauri star (Reipurth et al. 
1993; Chini et al. 1997). It has also been suggested that HH stars are more likely to be in 
binary systems (Reipurth 2000). The stellar and circumstellar properties inferred here may 
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be used to identify possible differences between stars which power HH flows and those which 
do not. 

The radial velocities of the HH stars newly measured here are consistent with the mean 
of Taurus, with the possible exception of HV Tau C (Section 3.1.1), which suggests that they 
are not more likely to be close (spectroscopic) binaries than non-HH stars. The locations 
of HH stars and non-HH stars on an H-R diagram are similar (Figure 14); HH stars do 
not appear to be systematically younger. We do find, however, a decreasing frequency of 
HH objects associated with the lowest mass stars and brown dwarfs (Sections 3.1.2 and 
3.2.2). The frequency may be even less, as none of the recently identified accreting T Tauri 
brown dwarfs (e.g. Muzerolle et al. 2003) are known to power HH flows. Although the 
strong forbidden line emission confirms that many low mass stars/BDs indeed power jets, 
the decreased frequency of spatially resolved regions of shocked emission suggests that their 
jets may be less powerful. This result is consistent with the mass dependent mass outflow 
rates and mechanical luminosities determined in Section 3.3.2. If the velocity imparted 
to the outflow also decreases with stellar mass (we have assumed a constant value in our 
calculations), this would lead to an even stronger mass dependence on the mechanical energy 
of the outflow. Although the vsmi values of HH stars are somewhat larger than those of 
non-HH stars (Section 3.1.1), this difference disappears if the values are normalized by their 
break-up velocity (Section 4.1). HH stars and non-HH stars have similar distirbutions of 
angular momentum. 

In comparison to non-HH stars, HH stars have slightly higher mass accretion rates (x2.5 
in the median), and much larger mass outflow rates (x20 in the median). It is unclear if 
the apparent difference in the mass outflow rate is a real time-averaged difference, or a 
consequence of recent accretion/outflow history. The similar line profiles of most HH stars 
and non-HH stars (Class I stars being the exception) suggests that edge-on orientation effects 
likely play less of a role in increasing the apparent forbidden line EWs (Section 3.1.4), which 
would translate into larger mass outflow rates. 

The distributions of millimeter fluxes, a direct tracer of circumstellar mass, for HH stars 
and non-HH stars are indistinguishable (Figure 21). This result is in contrast with the study 
of Reipurth et al. (1993) who found that HH energy sources generally have, on average, 
more than an order of magnitude 1.3 millimeter-wave flux densities than T Tauri stars in 
Taurus do. However, the Reipurth et al. (1993) comparison is biased in that the majority 
of HH stars in that study are much more luminous (1 — 10 3 L ) than the Taurus T Tauri 
stars compared to (~ 1 L Q ); as noted by Reipurth et al. (1993), the 1.3 millimeter flux 
shows a strong dependence on the bolometric luminosity. Our results suggest that when 
HH stars and non-HH stars of similar luminosities are compared, their millimeter fluxes 
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are likewise similar. The lower luminosities of Taurus HH stars may stem from current 
identification criteria (spatially resolved optical jet; Reipurth 1999). Less energetic jets and 
lower luminosity sources are easier to identify in Taurus given it's close proximity and low 
extinction. It is unclear if the more luminous HH stars in other star forming regions are at 
an earlier evolutionary stage, or simply more massive. 

Finally, we note that since HH objects are more often associated with Class I stars than 
Class II stars, the visual extinction of HH stars is systematically larger than for non-HH stars. 
Overall, however, with the exception of apparently larger mass outflow rates, the stellar and 
circumstellar properties of HH stars and non-HH stars of similar mass are similar. 

5. Summary 

Using the W. M. Keck I telescope, we obtained high dispersion (R ~ 34,000) optical 
(6330 - 8750 A) spectra of 15 Class I stars and 21 Class II stars in the nearby Taurus star 
forming region. Targets were selected based on evidence for either an infrared dominated 
luminosity as quantified by the SED diagnostics T\, \ and a, or the presence of a spatially 
resolved optical jet (Section 2.1). The optical emission from these environmentally young 
stars is spatially extended scattered light in some cases, and faint point-like in others. For 28 
of these 36 stars, our measurements are the first high dispersion optical spectra ever obtained. 
Photospheric features are detected in 11 of the Class I stars (42% of known Taurus Class 
I stars) and in all 21 of the Class II stars; strong emission lines (e.g. Ha) are detected in 
the spectra of all stars observed, even Class I stars not visible on the POSS-II red plates. 
Complementary J c -band images were obtained for the majority of Class I stars observed 
sp ectroscopically. 

Radial and rotational velocities are determined via cross-correlation with rotationally 
broadened spectral standards. All stars have radial velocities that are consistent (within 
3er) with the mean of Taurus. No spectroscopic binaries are identified, though HV Tau 
C and HK Tau B are noted as possible single-lined spectroscopic binaries. All stars are 
slowly rotating (fsim < 35 km/s). Spectral types and continuum excesses at 6500 A and 
8400 A are determined from the best fit rotationally broadened and veiled dwarf spectra. 
The inferred spectral types range from G8 to M6, with typical uncertainties of 1 spectral 
subclass. The new spectral types, in combination with J, if, K s photometry from the 2MASS 
database, are used to estimate visual extinctions and stellar luminosities. Masses and ages 
are determined from comparison with the Siess et al. (2000) evolutionary models. Surface 
gravity signatures in all spectra appear either dwarf-like or intermediate between those of 
dwarfs and giants, consistent with the modestly less than dwarf surface gravities expected 
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from their H-R diagram positions. Emission features associated with pressure sensitive lines 
inhibits direct determination of surface gravities from line profile analyses. 

The measured continuum excesses range from 0.0 to 5.7 times that of the photosphere, 
with one source (IRAS 04303+2240) showing variations by a factor of ~ 3 between 2 ob- 
servational epochs. For all stars, the continuum excesses appear to exhibit the trend of 
retaining constant flux over the 6500 - 8400 A wavelength interval. Although this has been 
observed previously (Basri & Batalha 1990), current hot spot models (e.g. Kenyon et al. 
1994; Calvet & Gullbring 1998) can not account for it. We interpret this emission as a 
cooler component of the shock, although emission from the inner disk can not be completely 
excluded. Mass accretion rates are determined from the excess emission at 6500 A under 
the assumption of a magnetically channeled accretion flow. Mass accretion rates for K7-M1 
spectral types span ~ 2 orders of magnitude, with a median value of 4 x 10~ 8 M Q /yr. This 
median value is larger than previous estimates (e.g. Gullbring et al. 1998) determined from 
excess emission at shorter wavelengths < 0.5 /im and based on a model that does not account 
for the observed excesses at red optical wavelengths. Until the continuum excess spectrum 
can be more accurately modeled, this will contribute a factor of 3 systematic uncertainty in 
mass accretion rates. 

Ha emission is detected from all sources, and in most CclS6S db wealth of permitted (Fe 
II, Ca II) and forbidden ([OI], [Nil], [SII]) emission lines are seen. Based on the relatively 
distinct emission line profiles of 3 optically veiled (and presumably accreting) edge-on disk 
systems (HH 30, HV Tau, HK Tau B), we suggest that CoKu Tau 1, IRAS 04260+2642, and 
ZZ Tau IRS also have nearly edge-on orientations. We confirm previous results that find 
larger forbidden-line emission associated with Class I stars than Class II stars. However, we 
attribute this to an orientation bias that allows a more direct view of the forbidden emission 
line region than the stellar photosphere, and not to larger mass outflow rates. Mass outflow 
rates are determined from the strength of [SII] 6731 A emission, under the assumption of 
a bipolar jet. Excluding edge-on disk systems, which may bias the measured forbidden-line 
EWs to larger values, the mass outflow rates for K7-M1 spectral types span ~ 3 orders of 
magnitude, with a median value of 2 x 10~ 9 M Q /yr. The ratios of M out fi ow /M in fi ow span 
~ 2 orders of magnitude, with a median value of 0.05. 

The inferred stellar and circumstellar properties are used to conduct statistical compar- 
isons of Class I and Class II stars, as well as HH stars and non-HH stars. A summary of 
the median statistics is compiled in Table 5, with brown dwarfs excluded. The distribution 
of stellar masses of Class I stars is similar to that of Class II stars, ranging from substellar 
to several solar masses. Of particular interest are IRAS 04158+2805, IRAS 04248+2612, 
and IRAS 04489+3042 which have substellar masses. These are the first spectroscopically 
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confirmed Class I brown dwarfs. Brown dwarfs as low in mass as 0.05 M Q likely form via 
dynamical collapse of a cloud core, experience an embedded Class I phase of evolution and 
are capable of powering molecular and HH flows. Stellar luminosities of optically revealed 
Class I stars suggest ages of ~ 10 6 yr, consistent with the ages of Class II and Class III T 
Tauri stars but inconsistent with ages implied by stellar birthline predictions and relative 
number statistics (~ 10 5 yr). In light of the wide range of disk dispersal time scales among 
T Tauri stars, we speculate that many Class I stars represent a T Tauri subsample with the 
longest envelope/disk dispersal time scale. 

For both Class I and Class II stars, approximately 25% of the bolometric luminosity 
is generated through disk accretion. The majority of the bolometric luminosity originates 
from the star. This result strongly supports one proposed resolution of the "luminosity 
problem" for Class I stars - they do not have accretion dominated luminosities. We further 
propose that most Class I stars in Taurus are past their main accretion phase. In some 
cases, Class I stars may actually be Class II stars that have been misclassified because of 
nearly edge-on orientations, or because of a biased SED caused by spatially unresolved 25 
/im measurements. In the more general case, however, we suggest that the envelope infall 
rates, which are roughly 2 orders of magnitude larger than the observed disk accretion rates, 
have been over-estimated by not properly accounting for the emission from a circumstellar 
disk. 

While the proposed scenario for Class I stars would explain the similar stellar proper- 
ties (masses, ages, rotation rates) and circumstellar properties (mass accretion rates, mass 
outflow rates) of Class I and Class II stars, there are still several issues that challenge this 
interpretation. The positions of Class I stars are correlated with the positions of dense cores 
(Hartmann et al. 2002) suggesting they are still close to their birth site. It is somewhat 
surprising that they are not more dispersed, like Class II stars, if Class I stars are indeed as 
old as Class II stars. The statistical ages of the few known Class stars suggest a very rapid 
formation timescale, and large mass accretion rates; many Class stars will have an even 
more severe luminosity problem (e.g. IRAS 04368+2557). Moreover, their rarity implies 
that there is very little ongoing star formation in Taurus. We emphasize, however, that the 
limited sensitivity of current surveys make it difficult to assess the validity of these so-called 
statistical problems. Observations with the Spitzer Space Telescope will likely yield a more 
comprehensive picture of star formation in the Taurus molecular cloud. 

Many of the Class I and Class II stars studied here have been identified as the energy 
sources of HH objects. The primary difference between HH stars and non-HH stars is that 
HH stars have stronger forbidden line emission, which translates into larger mass outflow 
rates by a factor of ~ 20 in the median. Whether this is a real time-averaged difference or a 
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consequence of recent accretion/outflow history is unclear. We also find that HH flows are 
less commonly associated with very low mass stars (< 0.2 M Q ) and brown dwarfs than with 
more massive stars. One remarkable exception is the Class 1 brown dwarf IRAS 04248+2612 
which powers both an HH and a molecular flow. Overall, with the exception of larger mass 
outflow rates, the stellar and circumstellar properties of HH stars and non-HH stars of similar 
mass are generally indistinguishable. 

On a concluding and reflective note, as we anticipate the exciting young star discoveries 
likely to be provided by the sensitive Spitzer Space Telescope, it is perhaps fitting that we are 
now determining fundamental stellar and circumstellar properties of environmentally young 
stars identified 20 years ago by Spitzer 's predecessor, the IRAS satellite. The extensive sur- 
veys of star forming regions being conducted with Spitzer may revolutionize our generally 
accepted ideas regarding star formation in the same way that the enlightening discoveries 
of the IRAS satellite did. However, proper interpretation of newly discovered red, faint, 
and potentially very young objects will likely depend critically on understanding the rela- 
tion between the stellar, accretion, and outflow properties, often best studied short-ward of 
~ 1/im, and the disk+envelope properties Spitzer will study at longer wavelengths. Obser- 
vations of stars in the Taurus star forming regions, with its close proximity and the low cloud 
extinction, will likely continue to play an important role in this endeavor. 

We are grateful to G. Doppmann, S. Edwards, J. Eisner, P. Hartigan, L. Hartmann, M. 
Liu and P. Williams for helpful discussions. This publication makes use of data products from 
the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts 
and the Infrared Processing and Analysis Center/ California Institute of Technology, funded 
by the National Aeronautics and Space Administration and the National Science Foundation. 
Finally, we recognize and acknowledge the very significant cultural role and reverence that 
the summit of Mauna Kea has always had within the indigenous Hawaiian community. We 
are most fortunate to have the opportunity to conduct observations from this mountain. 
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Fig. 1. — Spectral index versus bolometric temperature for stars in Taurus. A total of 44 
stars have either T bot < 650° Kora> 0.0 or power an HH flow. We consider these to be 
environmentally young stars (Table 1). The subset of Class I stars have both T boi < 650° 
K and a > 0.0 (when both values are available; see text). Large symbols indicate stars that 
power HH flows while small symbols represent stars that do not. Diamonds indicate stars 
that have been observed spectroscopically previously (Basri & Batalha 1990; Hartigan et al. 
1995; White & Basri 2003) and squares indicate stars observed in this study, 15 of which are 
within the Class I regime. Of the 36 stars observed here, 6 have no measurable continuum 
(unfilled squares; only 5 are visibly plotted since MHO-1 and MHO-2 overlap) while 30 (3 
of which are binary star components with overlapping positions) show a stellar continuum 
(filled squares) which we use to study the stellar and accretion properties. The X's indicate 
stars that have never been observed at high dispersion. 
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Fig. 2. — Spectral segments obtained with Keck/HIRES of environmentally young stars. 
Only stars with S/N > 2 are plotted; portions of spectra with S/N < 2 are shown as dashed 
lines. The 6445-6500 A, 7060-7130 A, 8420-8455 A regions are the primary ones used to 
determine the spectral types and veiling; the best fit dwarf spectral standards, rotationally 
broadened and veiled, are shown for each star. The youth diagnostic Li I 6708 A and the 
gravity sensitive feature KI 7699 A are also shown; these are typically the strongest and 
most easily identifiable lines in low S/N optically veiled spectra. Several optically veiled, 
accreting stars show Fe II 6456 A and OI 8446 A emission. 
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Fig. 3. — Same as Figure 2. 
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Fig. 4. — Spectral segments obtained with Keck/HIRES of additional T Tauri stars not 
meeting the environmentally young criteria. Features are as in Figure 2. Spectra from 2 
epochs are shown for IRAS 04303+2240 (see Table 2). This star was both brighter and more 
optically veiled during the first epoch. 




Fig. 5.— Ha, Ca II 8498A, [01] 6364 A, [Nil] 6583 A, [SII] 6731 A emission lines of the 
observed sample (Table 2). All spectra are shifted in wavelength to the mean systemic 
velocity of Taurus (18 km/s; vertical dashed line) and normalized to unity in the continuum. 
The vertical scale runs from 0.0 to 3.0, except when the peak to continuum ratio exceeds 
3.0, in which case the scale runs from 0.0 to the peak value. 
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Fig. 6. — Same as Figure 5. 
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Fig. 7. — Same as Figure 5. 
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Fig. 8. — Same as Figure 5. 
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Fig. 9. — Projected rotational velocities (vsiai) versus bolometric temperature (top panel) 
and spectral index (bottom panel). The vertical dotted lines are proposed values for separat- 
ing Class I and Class II stars (see Section 2.1). Large symbols indicate stars that power HH 
flows (ie. HH stars). Open symbols are i>sini upper limits. The stellar rotational velocities do 
not evolve significantly with either evolutionary diagnostic. HH stars rotate slightly faster 
than non-HH stars, on average. 
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Fig. 10. — Ratio of optical veiling at 8400 A to that at 6500 A versus stellar temperature. 
Squares are ratios determined from values listed in Table 2 and from White & Basri (2003); 
Diamonds are determined from the values in Hartigan & Kenyon (2003, r 6110 to r 8115 for 
these points). The observed veiling ratios are more consistent with a continuum excess of 
constant flux (F ex = C; dashed line) over this wavelength interval, than with a Rayleigh- 
Jeans approximation (F ex oc A -4 ; dotted line) as would be expected if the excess originated 
from a spot much hotter than the stellar photosphere. We interpret this as evidence for a 
second, cooler component to the accretion generated excess emission. 
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Fig. 11. — Optical veiling near 6500 A (r 6500 in Table 2) versus bolometric temperature (top 
panel) and spectral index (bottom panel). The symbols are the same as in Figure 9. The 
amount of excess emission does not evolve significantly with either evolutionary diagnostic. 
HH stars and non-HH stars have similar excesses. 
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Fig. 12. — Observed EW[Ha] versus bolometric temperature (top panel) and spectral index 
(bottom panel). The symbols are the same as in Figure 9; circled symbols indicate stars 
with known or likely edge-on disk systems (see Section 3.1.4). The EW[Ha] does not evolve 
significantly with either evolutionary diagnostic. HH stars and non-HH stars have similar 
EW[Ha] values. 
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Fig. 13. — Observed EW [SII 6731] emission versus bolometric temperature (top panel) and 
spectral index (bottom panel). The dashed lines and symbols are the same as in Figure 
9; circled symbols indicate stars with known or likely edge-on disk systems. The average 
EW[SII] decreases toward hotter bolometric temperatures and negative spectral indices. HH 
stars have systematically larger EW[SII]s than non-HH stars. Stars with edge-on orientation 
have EW[SII]s biased towards artificially large values (Section 3.1.4). 
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diagnostics: spectral index > 0.0, Tb i < 650, and stars that power HH flows. Values 
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the dispersion in these luminosities. The evolutionary models of Siess et al. (2000) are shown 
with isochrones at 10 6 , 10 7 , 10 8 years (solid lines), and mass tracks at 0.2, 0.5, 1.0, 2.0, 4.0 
M (dotted lines). The thick dot-dashed line indicates the "stellar birthline" for the case of 
Macc = 10 -5 M /yr (Fletcher & Stahler 1994). The environmentally young stars appear to 
have ages of ~ 10 6 years, similar to the mean ages of environmentally older Taurus T Tauri 
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Fig. 15. — Stellar mass versus bolometric temperature (top panel) and spectral index (bottom 
panel). The symbols are the same as in Figure 9. Class I stars have masses similar to 
Class II stars, ranging from substellar masses to several solar masses. Note that although 
many Class II brown dwarfs are known, they are not plotted here since they do not have 
bolometric temperature or spectral index values because of insufficient infrared and sub- 
millimeter measurements. 
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Fig. 16. — Projected rotational velocities (i>sini) normalized by the break-up velocity versus 
bolometric temperature (top panel) and spectral index (bottom panel). The symbols are 
the same as in Figure 9. Angular momentum does not evolve significantly with either 
evolutionary diagnostic. The angular momentum of HH stars and non-HH stars are similar. 
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Fig. 17. — Quantities used in the conversion of optical veiling values at 6500 A to mass 
accretion rates (veiling, stellar luminosity, excess flux in the 6000 - 6500 A passband, total 
accretion luminosity, mass accretion rate) versus stellar mass. Squares are new veiling mea- 
surements and diamonds are veiling measurements from the literature (see Section 3.3.1). 
The vertical error bar in each panel indicates the range of values determined by Gullbring 
et al. (1998) from a broader but bluer wavelength interval (3200 - 5300 A), constrained to 
the same mass range shown here. Surprisingly, the excess flux levels measured over the 6000 
- 6500 A passband are similar to those measured by Gullbring et al. (1998); the standard 
hot spot model predicts they should be only one-tenth. A factor of 2 can be attributed 
larger stellar luminosities, from which the excess fluxes are determined relative to, but the 
red excess fluxes are still a factor of ~ 5 greater than predicted. A reduced bolometric 
correction is adopted to convert the excess flux levels to total accretion luminosity, which 
leads to accretion luminosities onlv ~ 3 times larger than those of fGullbrine' et al. 19981. 
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Fig. 18. — Mass accretion rate versus bolometric temperature (top panel) and spectral index 
(bottom panel). The symbols are the same as in Figure 9. Although the mass accretion rates 
for Class I stars extend to lower values than for Class II stars, this is primarily a consequence 
of lower masses. Over a similar stellar mass range, the mass accretion rates do not evolve 
significantly with either evolutionary diagnostic. HH stars have mass accretion rates that 
are similar to those of non-HH stars, although larger by a factor of ~ 2.5 on average. 
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Fig. 19. — Mass outflow rate versus bolometric temperature (top panel) and spectral index 
(bottom panel). The symbols are the same as in Figure 9; circled symbols indicate stars with 
known or likely edge-on disk systems (see Section 3.1.4). Over a similar stellar mass range, 
the mass outflow rates appears to decline by a median factor of ~ 20 with either bolometric 
temperature or spectral index. Similarlly, HH stars have mass outflow rates that are, in the 
median, ~ 20 times that of non-HH stars. These difference may be a consequence of biases 
in the observed EW[SII] emission, as is the case for known edge-on disk systems 
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Fig. 20. — Ratio of mass outflow to mass accretion rate versus bolometric temperature (top 
panel) and spectral index (bottom panel). The symbols are as in Figure 9; upperward pointing 
triangles are lower limits and downward pointing triangles are upper limits. Circled symbols 
indicate stars with known or likely edge-on disk systems. The right hand ordinate shows 
the corresponding ratio of outflow mechanical lumionsity to accretion luminosity, assuming 
a 0.5 M Q star with a radius of 2.0 R . Class I stars and HH stars have systematically larger 
ratios than Class II stars or non-HH stars. These difference may be a consequence of biases 
in the observed EW[SII] emission, as is the case for known edge-on disk systems. 
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Fig. 21. — 1.3-mm Flux density versus bolometric temperature (top panel) and spectral index 
(bottom panel). The right hand ordinate shows the conversion to disk mass, as described 
in the text (Section 4.6). The symbols are as in Figure 9; the crosses indicate Class I 
stars without spectroscopically determined stellar properties. The disk mass does not evolve 
significantly with either evolutionary diagnostic. HH stars and non-HH stars have similar 
disk masses. 



Table 1. A Sample of Young Stars in Taurus- Auriga 
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+22 54 36.0 


(+0.38) 




I: 




yes 




11.14 


wh 


DO Tau 


04353+2604 


prev 


04 38 28.58 


+26 10 49.4 


(-0.51) 


(2125) 


II 


yes 


yes 


11.17 


9.47 


kh95 



Table 1 — Continued 



Common 


Associated 










Spectral 




Evol 


HH e 


visible^ 


Ic 


J 9 


jh 


Name 


IRAS Source 


Obs?" 


RA(2000) 6 


DEC(2000) 6 


Indcx c 


(°K) 


Class d 


source? 


in DSS? 


(mag) 


(mag) 


Rcl 


HV Tau C 


04353+2604 


new 


04 38 35.48 


+26 10 41.5 


(-0.51) 


(2125) 


II 


yes 


yes 


14.93 




mm94 




04361+2547 AB 




04 39 13.89 


+25 53 20.9 


(+1.54) 


(144) 


I 




no 


22.3: 


16.44 


wh 




04365+2535 




04 39 35.19 


+25 41 44.7 


+1.26 


172 


I 




no 


23.5: 


16.91 


wh 




04368+2557 




04 39 53.8 


+26 03 06 


+0.74 


59 





yes 


no 


19.0: 




wh 


1C 2087 


04369+2539 


new 


04 


39 55.75 


+25 45 02.0 


-0.57 


860 


II 


yes 


yes 




10.66 


wh 




04381+2540 




04 


41 12.68 


+25 46 35.4 


+1.29 


139 


I 


yes 


no 


21.1: 


17.15 


wh 


Haro 6-33 


04385+2550 


new 


04 


41 38.82 


+25 56 26.8 


+0.17 


636 


r 


yes 


yes 


15.58 


11.85 


kh95 


DP Tau 


04395+2509 


new 


04 


42 37.70 


+25 15 37.5 


-0.44 


1991 


ii 


yes 


yes 


11.95 


11.00 


kh95 


UY Aur AB 


04486+3042 


prev 


04 


51 47.38 


+30 47 13.5 


(-0.29) 


(2120) 


ii 


yes 


yes 


10.83 


9.13 


kh95 




04489+3042 


new 


04 


52 06.68 


+30 47 17.5 


+0.29 


399 


r 




no 


20.3: 


14.43 


wh 


RW Aur AB 


05046+3020 


new 


05 


07 49.54 


+30 24 05.1 


(-0.56) 


(4478) 


ii 


yes 


yes 


9.34 


8.38 


kh95 












Additional T Tauri Stars in 


Taurus-Auriga 














04278+2253 


new 


04 


30 50.28 


+23 00 08.8 


(-0.34) 


(856) 


ii 




yes 




8.78 


wh 


HK Tau 


04288+2417 


new 


04 


31 50.57 


+24 24 18.1 


(-0.51) 


(2148) 


ii 




yes 


12.37 


10.45 


kh95 


Haro 6-13 


04292+2422 


new 


04 


32 15.41 


+24 28 59.7 


-0.20 


910 


ii 




yes 




11.24 


wh 




04303+2240 


new 


04 


33 19.07 


+22 46 34.2 


-0.35 


886 


ii 




yes 


15.73 


11.10 


kh95 


DN Tau 


04324+2408 


new 


04 


35 27.37 


+24 14 58.9 


-0.84 


2890 


ii 


no 


yes 


10.49 


9.14 


kh95 


V836 Tau 




new 


05 


03 06.60 


+25 23 19.7 


-0.93 


3462 


ii 


no 


yes 


11.19 


9.91 


kh95 



indicates optical observations at high dispersion from this study ('new') or previous work ('prev'). 

b Positions are from the 2MASS Point Source Catalogue if available. The positions for 3 of the remaining 6 stars are from the 2MASS Extended Source Catalogue 
(IRAS 04016+2610, DC Tau B, and IRAS 04248+2612). The positions for IRAS 04166+2706 and IRAS 04368+2557, which arc too faint to be detected by the 
2MASS survey, are from Park & Kcnyon (2002) and Kcnyon ct al. (1990), respectively. The position of HV Tau C is determined from the 2MASS measurement 
of HV Tau AB offset by the separation and position angle measurements of Woitas & Leinert (1998). The position listed for the binary star MHO 1/2 is that of 
the primary, MHO 2. 

c Spectral indices are determined from the K-[25] colors listed in Kenyon & Hartmann (1995). Bolometric temperatures are from Chen et al. (1995). Values 
listed in parentheses have been determined from colors or spectral energy distribution that include contributions from more than 1 star. 

d The Class type (0, I or II) is set according to the bolometric temperature and spectral index criteria discussed in Section 2.1. Class I stars which are unresolved 
(ie. point-like) at 1.3 millimeters are marked with an apostrophe (Motte & Andre 2001); Haro 6-28, which is marked with a colon, has no spatial millimeter 
information. 

indicates if the star powers a Herbig-Haro flow (Gomez et al. 1997; Reipurth 1999). 
f Indicates if the object is visible in the Second Generation POSS-II Red plates (ref). 



sj magnitudes arc from the 2MASS Point Source Catalogue except for IRAS 04016+2610, DG Tau B, and IRAS 04248+2612, which are from the 1. 
Extended Source Catalogue. 

h kh95 = Kenyon & Hartmann (1995); wh = this work; mm94 = Magazzu & Martin (1994); mf83 = Mundt & Fried (1983) 
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Table 2. Observational Summary and Measured Photospheric Properties 



Exp. Radial 
Obs Time S/N S/N EW[LiI] Velocity usini 

Star Date (min) 6700 8450 (A) (km/s) (km/s) SpT r 650 o r 8 400 



Environmentally Young Stars in Taurus- Auriga 



IRAS 04016+2610 


99 


Dec 


5,6 


215 


1.6 


2.1 


< 0.5 


19 ± 7 : 


< 15 : 


K: 


< 0.3 : 


< 0.37 


MHO 1 


03 


Feb 


17 


10 


0.1 


1.3 


< 4.2 












MHO 2 


03 


Feb 


17 


30 


5.0 


18 


0.39 ± 0.06 


14.8 ± 3.0 


21.5 ± 4.2 


M3.5+1 


< 0.24 


< 0.05 


IRAS 04154+2823 


03 


Feb 


17 


60 


0.2 


0.9 


< 1.4 












CoKu Tau 1 


02 


Dec 


13 


30 


8.5 


15 


0.49 ± 0.03 


15.0 ± 0.8 


15.3 ± 1.7 


K7+1+M 


0.38 ± 0.18 


< 0.19 


IRAS 04158+2805 


02 


Dec 


13 


30 


4.0 


10 


0.21 ± 0.07 


18.1 ± 5.4 


23.1 ± 4.5 


M6+1 


< 1.0 


0.58 ± 0.30 


T Tau 


03 


Feb 


18 


1.5 


144 


119 


0.36 ± 0.02 


20.2 ± 0.4 


20.9 ± 1.0 


K0+2 


< 0.20 


< 0.19 


Haro 6-5 B 


02 


Dec 


13 


20 


7.0 


5.1 


0.45 ± 0.04 


18.6 ± 3.6 


20.7 ± 6.4 


K5+2 


1.0 ± 0.2 


1.4 ± 1.2 


DG Tau B 


02 


Dec 


13 


30 


0.6 


0.4 


< 0.5 












DG Tau 


99 


Dec 


6 


15 


59 


65 


0.24 ± 0.02 


19.3 ± 2.7 


28.6 ± 5.1 


K3+2 


2.1 ± 0.5 


1.8 ± 0.5 


IRAS 04248+2612 


99 


Dec 


6 


180 


1.3 


4.4 


< 0.2 


14.4 ± 2.5 


16.4 ± 4.2 


M5.5+1 


1.4 : 


0.60 ± 0.21 


IRAS 04260+2642 


02 


Dec 


13 


60 


5.3 


8.2 


0.43 ± 0.05 


22.3 ± 3.6 


< 18 


K6+2 


1.3 ± 0.5 


0.57 + 0.27 


GV Tau 


99 


Dec 


6 


30 


5.4 


7.1 


0.25 ± 0.05 


13.4 ± 3.8 


25.3 ± 7.0 


K7+2 


1.7 ± 0.4 


1.1 ± 0.4 


IRAS 04264+2433 


99 


Dec 


5 


60 


2.0 


5.1 


0.21 ± 0.14 


11 ± 8 : 


< 15 : 


Ml+2 


1.4 : 


0.99 ± 0.49 


ZZ Tau IRS 


03 


Feb 


18 


60 


4.7 


8.4 


0.15 + 0.06 


18.4 + 3.7 


22 ± 12 


M4.5+2 


2.1 + 0.7 


0.81 + 0.29 


L1551 IRS 5 


99 


Dec 


5 


60 


0.7 


0.1 


< 0.4 














03 


Feb 


17 


30 


0.6 


1.4 


< 0.5 












L1551 IRS 5 HH 


03 


Feb 


17 


30 


0.6 


0.1 


< 0.5 












L1551 IRS 5 Neb 


03 


Feb 


18 


60 


1.0 


0.3 


< 0.8 












HH 30 IRS 


02 


Dec 


13 


60 


9.6 


8.5 


0.17 + 0.03 


20.3 + 3.5 


< 12 


M0+2 


5.3 + 1.0 


5.7 + 0.9 


HL Tau 


99 


Dec 


5 


20 


29 


33 


0.26 + 0.02 


18.3 ± 2.6 


26.3 + 5.2 


K5+1 


0.91 ± 0.07 


0.99 + 0.15 


IRAS 04295+2251 


99 


Dec 


5 


60 


0.1 


0.8 


< 2.7 












IRAS 04302+2247 


99 


Dec 


5 


120 


0.8 


1.8 


< 0.4 












Haro 6-28 


02 


Dec 


13 


20 


28 


47 


0.46 ± 0.02 


16.1 + 0.3 


10.1 ± 1.0 


M3. 5+0.5 


0.23 ± 0.05 


< 0.07 


HV Tau C 


02 


Dec 


13 


40 


21 


19 


0.47 + 0.02 


23.1 + 0.6 


19.4 + 1.9 


K6+1 


0.39 ± 0.07 


< 0.19 


IC 2087 


02 


Dec 


13 


20 


2.5 


14 


< 0.12 


22 + 8 : 


< 15 : 


K: 


2.7 : 


3.2 + 0.6 


Haro 6-33 


02 


Dec 


13 


20 


14 


23 


0.59 + 0.02 


15.8 + 0.6 


22.8 + 1.6 


M0.5+0.5 


0.13 ± 0.06 


< 0.04 


DP Tau 


03 


Feb 


17 


5 


47 


48 


0.28 + 0.02 


16.8 + 0.6 


19.2 ± 1.3 


M0+1 


1.8 + 0.2 


1.3 + 0.3 


IRAS 04489+3042 


99 


Dec 


5 


60 


0.5 


2.5 


< 0.6 


15 + 5 : 


< 15 : 


M6+2 


1.3 : 


0.56 + 0.20 


RW Aur A 


03 


Feb 


17 


1 


67 


49 


0.29 + 0.02 


16.0 + 1.7 


34.4 + 5.7 


K2+2 


1.1 + 0.3 


1.5 + 0.6 


RW Aur B 


03 


Feb 


17 


1 


14 


16 


0.50 ± 0.02 


15.9 + 0.5 


12.2 ± 1.6 


K6+1 


0.26 ± 0.12 


0.17 + 0.13 



Additional T Tauri Stars in Taurus-Auriga 



Table 2 — Continued 



Exp. Radial 







Obs 




Time 


S/N 


S/N 


EW[Li I] 


Velocity 


usim 








Star 


Date 




(min) 


6700 


8450 


(A) 


(km/s) 


(km/s) 


SpT 


^6500 


^8400 


IRAS 04278+2253b 


02 


Dec 


13 


15 


14 


24 


0.42 ±0.02 


15.7 ±0.4 


< 7.4 


K7/M0±l 


0.61 ±0.10 


0.49 ±0.11 


IRAS 04278+2253a 


02 


Dec 


13 


15 


15 


27 


0.25 ±0.02 


19.2 ±2.7 


15.7 ±4.9 


G8±2 


0.34 ±0.22 


0.63 ±0.33 


HK Tau A 


02 


Dec 


13 


8 


23 


31 


0.50 ±0.02 


21.6 ±0.5 


17.8 ± 1.1 


Ml±0.5 


0.40 ±0.10 


0.24 ±0.06 


HK Tau B 


02 


Dec 


13 


20 


12 


16 


0.58 ±0.02 


17.4 ±0.3 


< 7.3 


Ml±0.5 


0.14 ±0.06 


< 0.12 


Haro 6-13 


02 


Dec 


13 


20 


24 


40 


0.50 ±0.02 


18.6 ±0.8 


23.3 ± 1.1 


M0±0.5 


0.43 ±0.08 


0.26 ±0.08 


IRAS 04303+2240 


99 


Dec 


6 


30 


53 


69 


< 0.05 








14.9 : 


7.9 : 




03 


Feb 


17 


20 


10.3 


18 


0.25 ±0.03 


20.8 ±3.6 


35 ± 16 


M0.5±l 


3.9 ±0.4 


2.6 ± 1.1 


DN Tau 


99 


Dec 


6 


20 


49 


55 


0.53 ±0.02 


16.9 ±0.3 


9.8 ±0.9 


M0±0.5 


0.24 ±0.07 


< 0.07 


V836 Tau 


99 


Dec 


6 


20 


71 


71 


0.55 ±0.02 


18.5 ±0.3 


12.1 ± 1.0 


K7±1±M 


0.12 ±0.10 


< 0.16 












T Tauri stars 


in the TW Hydrae Association 








TWA 8a 


03 


Feb 


18 


20 


195 


231 


0. ± 0.02 


8.7 ± 0.2 


< 5.5 


M3±0.5 


< 0.08 


< 0.08 


TWA 8b 


03 


Feb 


18 


20 


44 


88 


0. ± 0.02 


8.2 ±0.2 


10.8 ± 1.0 


M5.5±0.5 


< 0.06 


< 0.17 


TWA 9a 


03 


Feb 


18 


18 


239 


190 


0. ± 0.02 


10.9 ±0.3 


11.1 ± 1.0 


K7±l 


< 0.09 


< 0.10 


TWA 9b 


03 


Feb 


18 


18 


63 


88 


0. ± 0.02 


12.2 ±0.4 


9.0 ± 1.0 


M3.5±0.5 


< 0.20 


< 0.15 
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Table 3. Summary of Median Properties, Substellar Objects Excluded 

















Ha 10% 


















■usini 








EW[Ha] 




width 




EW[SII] 




\ogM A cc 




\ogM Qut 




type 


(km/s) 


N 


»"6500 


N 


(A) 


N 


(km/s) 


N 


(A) 


N 


(M /yr) 


N 


(M /yr) 


N 


Class Is 


19.8±5.7 


8 


0.6±0.6 


8 


-51±33 


13 


344 ± 131 


13 


-2.4±5.1 


13 


-7.1 ±0.8 


8 


-7.0 ±1.5 


8 


Class lis 


13.7±11.0 


41 


0.6±1.3 


50 


-59±48 


49 


429 ± 108 


46 


-0.1±17 


42 


-7.3 ±0.7 


51 


-8.2 ± 1.1 


42 


HH 


19.4±6.0 


19 


1.0±1.4 


23 


-61±69 


24 


409 ± 120 


21 


-1.8±22 


21 


-7.0 ±0.6 


23 


-7.2 ±0.9 


20 


non-HH 


11.4±11.9 


31 


0.5±1.1 


34 


-58±37 


37 


414 ± 116 


38 


-0.1±3.4 


33 


-7.4 ±0.7 


35 


-8.7 ±1.0 


30 
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Table 4. Permitted and Forbidden Emission-Line Properties 





EWHtt 


Ha 


EWCall 


EW Call 


EW [O I] 


EW [N II] 


EW [S II] 


EW [S II] 




6563 


10%- width 


8498 


8662 


6364 


6583 


6716 


6731 


Star 


(A) 


(km/s) 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 






Environmentally Young Stars in Taurus-Auriga 








TT) AO A ,d A 1 f i c\r--\ r\ 

IKAb U4Ulo+^olU 


41 


550 


4.4 


1.9 


0.6 


2.1 


l.i 


1.8 


T\ , T TU r\ 1 


58 


253 


17 


15 


< 4. 


< 4. 


< 4. 


, A 

< 4. 


TV ,T TT f~\ O 


59 


249 


0.61 


0.46 


1.0 


0.91 


0.7 


1.2 


T T") AO fl<1 r/ I O O O O 

IKAb U4154+z82o 


> 17 


340 






< 1.4 


< 1.4 


< 1.4 


< 1.4 


/~1 „ TV" , , rp -i 

Oolvu lau 1 


70 


285 


0.37 


0.50 


16 


19 


12 


21 


TTD AO n/i 1 i oonc; 
IKAb U41oo+z5Uo 


1/5 


287 


9.3 


6.8 


13 


6.1 


4.0 


8.6 


rp rp 

1 lau 


60 


380 


12 


12 


> 0.1 


0.11 


> 0.05 


0.41 


TT _ _ a r TD 

Haro d-o H 


91 


481 


31 


32 


< 0.1 


0.37 


< 0.04 


< 0.04 


DG Tau B 


> 276 


330 


> 10 


> 12 


> 36 


> 96 


> 298 


> 426 


p\ f~i rp 

JJLr lau 


63 


517 


49 


41 


3.0 


0.28 


1.1 


2.4 


TTD AO r\/iojo i o / ■ i 

IRAS 04248+2612 


163 


410 


14 


11 


8.5 


4.5 


6.8 


11 


ttd A o r\ 1 o ri i o£? i o 


126 


139 


4.0 


3.7 


19 


13 


8.9 


19 


GV Tau 


86 


552 


15 


13 


13 


3.1 


5.3 


8.0 


TTD AO rHO^yl 1 OJOO 

IKAb 04264+2433 


96 


367 


8.3 


7.5 


8.9 


4.9 


6.4 


9.8 


r/r/ rp TTD G 

ZZ lau 1Kb 


238 


268 


23 


19 


51 


30 


40 


77 


t irri TTj o c 

L1551 1Kb 5 


> 999 


196 


> 32 


> 44 


> 611 


> 999 


> 669 


> 935 


T i rri TTD O C 

L1551 1Kb 5 


> 412 


290 


< 0.8 


< 0.9 


> 281 


> 295 


> 127 


> 225 


T 1 KK1 TTD O C TTTJT 

L1551 1Kb 5 HH 


> 818 


196 


> 15 


> 18 


> 163 


> 685 


> 530 


> 739 


T i rri TTD O C XT^l-, 

L1551 1Kb 5 INcb 


44 


275 






< 1.7 


< 1.7 


< 1.7 


< 1.7 


TTTT on TTD O 

HH 1Kb 


199 


293 


13 


11 


49 


27 


58 


75 


HL Tau 


43 


390 


26 


22 


2.4 


1.2 


1.6 


2.4 


TTD AO f\ A O t \ f i o o r* 1 

1HAS 04295+2251 


> 11 


160 


> 3.3 


> 2.8 


< 3. 


< 3. 


< 3. 


< 3. 


TTZ> A Q H/IQnO 1 00/1*7 

IKAb 04o02+224/ 


6 / 


485 


36 


43 


< 2. 


< 3. 


2.8 


3.5 


TT ("JO 

Haro 0-28 


51 


344 


0.35 


0.24 


0.25 


< 0.1 


< 0.04 


< 0.04 


HV Tau C 


52 


261 


1.5 


1.0 


7.3 


10.3 


7.3 


13 


IC 2087 


22 


610 


14 


9.0 


< 0.4 


< 0.3 


< 0.1 


< 0.1 


Haro 6-33 


15 


377 


0.48 


0.31 


0.62 


1.0 


> 0.6 


1.5 


DP Tau 


74 


601 


13 


10 


1.8 


0.8 


0.85 


1.9 


IRAS 04489+3042 


> 96 


415 


11 


7.9 


> 7.2 


> 7.0 


> 2.7 


> 4.7 


RW Aur A 


42 


365 


36 


29 


< 0.13 


< 0.1 


> 0.03 


0.10 


RW Aur B 


17 


466 


5.9 


3.5 


0.11 


0.35 


0.16 


0.28 






Additional T Tauri Stars in Taurus- Auriga 








IRAS 04278+2253 B 


54 


498 


17 


15 


0.18 


0.12 


< 0.04 


< 0.04 


IRAS 04278+2253 A 


14 


434 


4.9 


3.8 


0.21 


0.10 


0.06 


0.14 


HK Tau A 


34 


435 


0.36 


0.42 


< 0.1 


< 0.1 


< 0.04 


< 0.04 


HK Tau B 


9.7 


196 


0.68 


0.59 


> 0.6 


0.53 


1.0 


1.7 


Haro 6-13 


34 


321 


1.8 


1.2 


1.1 


0.18 


> 0.15 


0.58 


IRAS 04303+2240 


67 


487 


30 


26 


0.35 


0.05 


0.14 


0.20 


IRAS 04303+2240 


137 


337 


24 


22 


9.5 


1.9 


2.4 


5.5 


DN Tau 


18 


291 


0.64 


0.54 


< 0.4 


< 0.4 


< 0.4 


< 0.4 


V836 Tau 


25 


270 


0.69 


0.47 


0.18 


< 0.4 


< 0.4 


< 0.4 
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Table 5. Inferred Stellar, Accretion and Outflow Properties 





Adopt 




A v 








1 Myr Mass 


log(MAcc) 


log(M 0llt ) 


log 


Star 


SpT 




(mag) 


(L©) 


(L©) 


(L©) 


(M ) 


(Mo /yr) a 


(Mo/yr) 


(Mn„t/MArr) 

\ LLL t /ILL / 










Environmentally Young Stars in 


Taurus-Auriga 








IRAS 04016+2610 


K4 


4580 


10.2 


0.45 


4.63 


3.70 


1.61 


-7.15 


-6.93 


+0.22 


MHO 1 












(> 1.60) 










MHO 2 


M3.5 


3180 


18.1 


2.61 


0.44 


(> 1.60) 


0.24 


< -8.48 


-8.99 


> -0.51 


IRAS 04154+2823 












> 0.25 










CoKu Tau 1 


K7 


4000 


6.80 


0.15 


1.60 


> 0.29 


0.68 


-7.36 


-6.41: 


+0.95: 


IRAS 04158+2805 


M6 


2760 


8.63 


0.050 


0.022 


> 0.44 


0.05 


< -9.50 


-9.72 


> -0.22 


T Tau 


K0 


5235 


4.24 


19.9 


29.6 


15.5 


3.38 


< -6.74 


-6.96 


> -0.21 


Haro 6-5 B 


K5 


4395 


9.96 


0.047 


3.0 


> 0.02 


1.20 


-6.76 


< -8.60 


< -1.84 


DG Tau B 












> 0.02 










DG Tau 


K3 


4775 


3.32 


3.62 


7.7 


6.36 


2.22 


-6.13 


-6.19 


-0.06 


IRAS 04248+2612 


M5.5 


2845 


7.02 


0.27 


0.044 


0.36 


0.07 


-8.97 


-9.11 


-0.13 


IRAS 04260+2642 


K6 


4200 


9.24 


0.054 


1.95 


0.09 


0.89 


-6.79 


-6.09: 


+0.70: 


GV Tau 


K7 


4000 


12.1 


1.82 


1.60 


6.98 


0.68 


-6.71 


-6.54 


+0.17 


IRAS 04264+2433 


Ml 


3605 


10.4 


0.14 


0.80 


0.37 


0.41 


-7.11 


-6.98 


+0.13 


ZZ Tau IRS 


M4.5 


3015 


7.60 


0.13 


0.21 




0.16 


-8.06 


-7.35: 


+0.70: 


L1551 IRS 5 












> 21.9 










HH 30 IRS 


M0 


3800 


2.96 


0.0064 


0.98 




0.52 


-6.45 


-5.48: 


+0.97: 


HL Tau 


K5 


4395 


7.43 


1.53 


3.0 


6.60 


1.20 


-6.80 


-6.84 


-0.04 


IRAS 04295+2251 












0.44 










IRAS 04302+2247 












0.34 










Haro 6-28 


M3.5 


3175 


4.38 


0.31 


0.44 


0.96 


0.24 


-8.50 


< -10.47 


< -1.97 


HV Tau C 


K6 


4200 






1.95 




0.89 


-7.31 


-6.46: 


+0.85: 


IC 2087 


K4 


4580 


18.1 


20.7 


4.6 


3.80 


1.61 


-6.20 


< -7.65 


< -1.46 


Haro 6-33 


M0.5 


3700 


10.2 


0.76 


0.87 


> 0.36 


0.46 


-8.11 


-8.03 


+0.08 


DP Tau 


M0 


3800 


6.31 


0.68 


0.98 


0.70 


0.52 


-6.92 


-7.42 


-0.51 


IRAS 04489+3042 


M6 


2760 


17.5 


0.24 


0.022 


0.30 


0.05 


-9.39 


> -9.92 


< -0.53 


RW Aur A 


K2 


4955 


1.58 


3.37 


12.9 


(3.20) 


2.84 


-6.32 


-7.63 


-1.31 


RW Aur B 


K6 


4200 






3.0 


(3.20) 


1.20 


-7.34 


-7.95 


-0.61 










Additional T Tauri Stars in Taurus-Auriga 








IRAS 04278+2253 B 


K8 


3900 






1.14 


(> 4.50) 


0.59 


-7.32 


< -9.24 


< -1.92 


IRAS 04278+2253 A 


G8 


5445 


11.4 


29.6 


41.9 


(> 4.50) 


3.54 


-5.18 


-7.08 


-0.90 


HK Tau A 


Ml 


3605 


5.42 


0.84 


0.80 


0.81 


0.41 


-7.65 


< -9.60 


-1.95 


HK Tau B 


Ml 


3605 






0.80 


> 0.01 


0.41 


-8.11 


-8.06: 


+0.05: 


Haro 6-13 


M0 


3800 


11.9 


2.11 


0.98 


1.30 


0.52 


-7.54 


-8.24 


-0.70 


IRAS 04303+2240 


M0.5 


3700 


11.7 


2.20 


0.87 


> 2.60 


0.46 


-6.05 


-7.76 


-1.71 


IRAS 04303+2240 


M0.5 


3700 


11.7 


2.20 


0.87 


> 2.60 


0.46 


-6.63 


-6.83 


-0.20 


DN Tau 


M0 


3800 


1.89 


0.68 


0.98 


1.00 


0.52 


-7.79 


< -9.46 


< -1.67 


V836 Tau 


K7 


4000 


1.68 


1.21 


1.60 


0.51 


0.68 


-7.86 


< -9.22 


< -1.36 



a Values marked with a colon indicate edge-on disk systems. The [SB] emission in these systems appears to be biased towards large values, lea< 
to overestimates of Mput an <J log(Mout/ Ma C c)- 



